EPJ manuscript No. 

(will be inserted by the editor) 



Small x Phenomenology - Summary and status 

The Small x Collaboration 

C*") Jeppe Andersen 1 , Serguei Baranov 2 , John Collins 3 , Yuri Dokshitzer 4 , Lidia Goerlich 5 , Giinter Grindhammer 6 , 

Gosta Gustafson 7 , Leif Jonsson 8 , Hannes Jung 8 , Jan Kwiecihski 5 ^ , Eugene Levin 9 , Artem Lipatov 10 , Leif Lonnblad 7 , 
y Misha Lublinsky 9 , Martin Maul 7 , Izabela Milcewicz 5 , Gabriela Miu 7 , Grazyna Nowak 5 , Torbjorn Sjostrand 7 , 

, Anna Stasto 11 , Nicusor Timneanu 12 , Jacek Turnau 5 , and Nikolai Zotov 13 
O . 

D ■ DAMTP and Cavendish Laboratory, University of Cambridge, UK, 



M. I 



Penn State Univ., 104 Davey Lab., University Park PA 16802, USA 

LPTHE, Universites P. & M. Curie (Paris VI) et Denis Diderot (Paris VII), Paris, France 
H. Niewodniczanski Institute of Nuclear Physics, Cracow, Poland 
Max Planck Institut, Munich, FRC, 

Department of Theoretical Physics, Lund University, Sweden 



' 8 Department of Physics, Lund University, Sweden 

, 9 School of Physics, Tel Aviv University, Tel Aviv, Israel, 

10 Moscow State University, Moscow, Russia 

' 11 DESY, Hamburg, FRC and H. Niewodniczanski Institute of Nuclear Physics, Cracow, Poland, 

t-H , 12 University of Uppsala, Sweden 

(f) i 13 Skobeltsyn Institute of Nuclear Physics, Moscow State University, Moscow, Russia 



' deceased 

23. Dec. 2003 



Oh! 



Abstract. A second workshop on small x physics, within the Small x Collaboration, was held in Lund 
in June 2002 with the aim of over-viewing recent theoretical progress in this area and summarizing the 
experimental status. 

This paper is dedicated to the memory of Jan Kwieciriski, who died unexpectedly on August 29, 2003. 
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1 Introduction as well as for detailed comparison with experimental data 

and the usage in the experimental groups at HERA and 
This paper is a summary of the 2nd workshop on small- a; elsewhere. Since then they have been used in very differ- 
parton dynamics held in Lund in the beginning of June ent areas, like jet production and heavy flavor physics. The 
2002. During two days we went through a number of theo- small- a; improved unintegrated parton densities obtained 
retical, phenomenological as well as experimental aspects from CCFM evolution implemented in the Monte Carlo 
of small- a; physics in short talks and long discussions, generators have been proven to be a very powerful tool 
Whereas our first workshop in 2001 and the resulting sum- in describing experimental data as well as for estimating 
mary paper [1] was dedicated to a general survey and dis- the effect of higher order corrections. For example only by 
cussion of small- a; physics in order to identify the most also applying the Cascade Monte Carlo in the extrac- 
pending questions, we concentrate here on those aspects, tion of F% and in the calculation of bottom production at 
where progress has been made, as well as on a more de- HERA, it was recognized that the extrapolation from the 
tailed discussion and some aspects of the experimental measured visible range to the total cross section is danger- 
situation. For a general introduction to small- a: physics ous and introduces large model dependencies. Now, in the 
and the small- a; evolution equations, as well as tools for area of bottom production, the visible cross sections are 
calculation in terms of Monte Carlo programs, we refer in reasonable agreement both with calculations applying 
the reader to [1]. fcx-factorization with CCFM evolved unintegrated gluon 
With the successful completion of the two full hadron density as well as with NLO calculations in the collinear 
level Monte Carlo programs LDCMC [2-5] and Cas- approach. This shows the importance of applying alterna- 
CADE [6,7], the necessary tools were provided for detailed 
studies both on a theoretical and phenomenological level 
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tive approaches even when extracting experimental mea- 
surements. 

This paper is organized as follows: First we discuss 
in more detail the definition of unintegrated gluon densi- 
ties, as well as the question on gauge invariance of parton 
densities in general and especially of the fcj_-factorization 
approach. In the following section we discuss results and 
problems in theoretical applications of the unintegrated 
parton distribution functions, different parameterizations, 
the scale in a s , the role of the non-singular terms in the 
g — > gg splitting function, saturation and the effects of en- 
ergy momentum conservation in the BFKL equation. The 
section also contains a discussion of polarized unintegrated 
distributions, polarization effects and color octet contribu- 
tions in J/tp meson production. The second part of this 
paper deals with experimental investigations of small- x 
effects and with the question, whether and where devia- 
tions from the collinear approach can be established, and 
whether a sign for a new evolution scheme like BFKL/ 
CCFM/ LDC has already been seen. We end this paper 
with an outlook and a definition of the next steps and 
goals. 



2 k± - factorization formalism 

The DGLAP[8-11] evolution treats successive parton 
emissions which are strongly ordered in virtuality and re- 
sums the resulting large logarithms of ratios of subsequent 
virtualities. Because of the strong ordering of virtualities, 
the virtuality of the parton entering the hard scattering 
matrix element can be neglected (treated collinear with 
the incoming hadron) compared to the large scale Q 2 . 

At very high energies, it is believed that the theoret- 
ically correct description is given by the BFKL [12-14] 
evolution. Here, each emitted gluon is assumed to take a 
large fraction, 1 — z\ z ^ of the energy of the propagating 
gluon, and large logarithms of 1/z are summed up to all 
orders. 

The CCFM [15-18] evolution equation resums also 
large logarithms of 1/(1 — z) in addition to the 1/z ones. 
Furthermore it introduces angular ordering of emissions 
to correctly treat gluon coherence effects. In the limit of 
asymptotic energies, it is almost equivalent to BFKL [19- 
21], but also similar to the DGLAP evolution for large x 
and high Q 2 . The cross section is fc^-factorized into an 
off-shell matrix element convoluted with an unintegrated 
parton density (uPDF) , which now also contains a depen- 
dence on the maximum angle 5 allowed in emissions. This 
maximum allowed angle S is defined by the hard scatter- 
ing quark box, producing the (heavy) quark pair and also 
defines the scale for which parton emissions are factorized 
into the uPDF. 

The original CCFM splitting function is given by: 

P ( z . a 2 k 2 )- ^(g'd-^) 2 ) , «s(fcL) A , -2 ,2 x 

1 Z{ Zi 

(1) 



with a s = and the non-Sudakov form factor A ns given 
by: 

, A f -2 . 2 x f 1 dz' f dq 2 _ 

In A ns ( Zi , qi , k ±i ) = - I —J-^a s - 

0{k ±l - q)G{q - z'qi) (2) 

fUz' dq 2 _ 
= - — — «s 3 

J Zi z J {z'qi)* T 

The angular ordering condition is given by: 

Zi-m-i < q~i (4) 

where the rescaled transverse momenta % of the emitted 
gluons is defined by: 



Here Zi — Xi/xi-i is the ratio of the energy fractions in the 
branching (i— 1) — > i and p±i is the transverse momentum 
of the emitted gluon i. The transverse momentum of the 
propagating gluon is given by It is interesting to note, 
that the angular ordering constraint, as given by eq.(4), 
reduces to ordering in transverse momenta p± for large 
z, whereas for z — > 0, the transverse momenta are free to 
perform a so-called random walk. 

In [1] it has been proposed to include also the non- 
singular terms in the splitting function as well as to 
consistently use fi r = p± for the renormalization scale 
in a s (fi r ), everywhere. These changes, although formally 
sub-leading, have significant influence for calculation per- 
formed at present collider energies. 

The inclusion of non-singular terms, as well as the evo- 
lution of quarks, is straightforward in the LDC model [2- 
5], which is a reformulation of CCFM, where the separa- 
tion between the initial- and final-state emissions is rede- 
fined. In addition to the angular ordering in eq.(4), the 
gluons emitted in the initial-state are required to have 

p±i > min(fc_u,fc_u_i). (6) 

In the double leading logarithmic approximation (DLLA) , 
this requires a reweighting of each splitting, completely 
canceling the non-Sudakov form factor, reducing the split- 
ting function in (1) to the leading singularities of the stan- 
dard DGLAP one, making the inclusion of non-singular 
terms as well as quark splittings a trivial exercise. The 
constraint in (6) means that the p± of the emitted gluon 
is always close to the highest scale in the splitting and the 
argument in a s is naturally taken to be p\ . 

While formally equivalent to the DLLA accuracy for 
the inclusive observable F 2 , it is important to note that 
the sets of chains of initial-state splittings summed over, 
arc different in LDC and CCFM. Therefore results for 
exclusive final states agree only after addition of final state 
radiation in the appropriate kinematical regions (which 
are different in the two formalisms). 

We here also want to mention the formalisms devel- 
oped in Rcfs. [22,23]. An evolution equation for a sin- 
gle scale uPDF, which interpolates between DGLAP and 
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BFKL, is presented by Kwiecinski, Martin and Stasto in 
[22]. The formalism for a two-scale uPDF by Kimber, Mar- 
tin and Ryskin [23] is based on the same single scale evo- 
lution equation, but an angular cut is applied for the last 
step in the chain. 



2.1 Unintegrated parton distributions 

In the following we discuss in detail the precise definition 
of (integrated or unintegrated) parton density functions 
(PDFs) [24]. 

1. A PDF is not a physical quantity in and of itself. It 
merely is a useful tool. 

2. A PDF is often given as a probability density of quarks 
or gluons within the framework of light-front quantiza- 
tion of QCD in the light-cone gauge. Such a definition 
is useful to provide motivation, intuition and an ini- 
tial candidate for a formal definition. But this method 
does not necessarily provide a valid definition. 

3. Whether or not some kind of consistent probability 
interpretation can be made with modified definitions 
is an open question. For many applications of PDFs, 
the answer to this question is irrelevant. 

4. The physical significance of PDFs is that there are 
factorization formulae involving them 1 . Factorization 
formulae (in their most general sense) give approxima- 
tions to physical amplitudes or cross sections that are 
useful and predictive because: 

(a) The PDFs are universal - the same in a range of 
different processes and situations. 

(b) Some (not necessarily all) of the coefficients in a 
factorization formula may be estimated, for exam- 
ple in fixed-order perturbation theory with a weak 
coupling. 

(c) Kernels of evolution equations (DGLAP etc) may 
similarly be estimated perturbatively. 

5. Since a PDF will include non-perturbative physics, it is 
generally desirable that an explicit definition be given, 
for example in terms of some Green function or a ma- 
trix clement of some (usually non-local) operator. 

6. Given point 1, it is not necessary that a PDF's defini- 
tion is explicitly gauge invariant. However, if the defini- 
tion is not gauge-invariant, the choice of gauge must be 
explicitly specified. It should be possible to transform 
the PDF's definition into an explicitly gauge-invariant 
form. But in general there should be extra parame- 
ters) for the parton density corresponding perhaps to 
a gauge-fixing vector or the direction of Wilson line 
factors in the operators. It will also be necessary to 
obtain evolution equations with respect to the extra 
variable(s). See the work of Dokshitzer, Diakonov and 
Troian [25] and of Collins and Soper [26,27] for exam- 
ple. 



7. The most obvious candidate definition for a PDF is as 
a number density in the light-cone gauge, essentially 



f(x,k±) 



(P\bjb k \p) 
(p\p) 



(7) 



where b\ and b k are creation and annihilation opera- 
tors for a flavor of parton in the sense of light-front 
quantization in light-cone gauge. However, as we will 
see below, such a definition is divergent beyond the 
lowest order of perturbation theory. 

8. The divergence arises from an integral over rapidity of 
emitted gluons and is present even if all IR and UV di- 
vergences are regulated. The divergence is an endpoint 
divergence due to the l/k + singularity in the gluon 
propagator in light-cone gauge: J dk + /k + Therefore 
it cannot be removed by a modification of the inte- 
gration path, and in that way changing the analytic 
prescription of the singularity. 

9. For an uPDF, the divergence cannot be canceled be- 
tween real final state and virtual gluon emission: Vir- 
tual gluon emission has an unrestricted transverse mo- 
mentum integral, but real gluon emission is restricted 
by the transverse momentum of the emitting parton 
(Fig. 1). Hence a cancellation of real and virtual diver- 
gences cannot occur simultaneously for all values of the 
transverse momentum of the emitted parton. This is a 
problem because without a cancellation between real 
and virtual divergencies the resulting parton density 
function diverges and becomes meaningless. 



virtual 
correction 



real emission 
P 2 




1 The status of a given factorization formula may be any- 
where from being a completely proved theorem to merely being 
a conjecture. 



Fig. 1. A schematic drawing indicating the real and virtual 
corrections in a parton splitting process. The transverse mo- 
mentum of the emitted parton p2 is limited by the emitting 
parton pi , whereas the transverse momenta in the virtual loops 
can reach any value. 



10. The rapidity divergence (J dk+/k + = J_ dy) involves 
gluon momenta in a region that has no relevance to 
the process: The momenta have infinite rapidity rel- 
ative to physical particles. Any sensible definition of 
a PDF must have a cutoff. A simple candidate would 
be obtained by taking the (incorrect) light-cone gauge 
definition but with the use of a planar gauge: n- A = 0, 
with n 2 7^ 0. The unintegrated parton density then 
has two extra parameters beyond the x and k± kine- 
matic variables. These are the inevitable renormaliza- 
tion scale /i and the variable (p ■ n) 2 /n 2 p 2 . The renor- 
malization scale fi has the approximate interpretation 
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of a cut off on the transverse momentum or the virtu- 
ally of virtual particles, and the last variable equals 
cosh 2 y, where y is the rapidity difference between the 
target and the gauge fixing vector. 

f f . In the CCFM formalism these two extra parameters 
are correlated. Thus the CCFM parameter q deter- 
mines both the transverse momentum cutoff and the 
limiting rapidity in the p-rest frame through the rela- 
tion y = ln(q/m p Xg). 

12. In Rcf. [24] Collins has proposed explicit gauge- 
invariant definitions of unintegrated PDFs that avoid 
the difficulties mentioned above. The evolution equa- 
tion would be that given by Collins, Soper and Ster- 
man [26,27]. Details in this approach are being worked 
out. 



2.2 Further questions on gauge invariance 

The question of gauge invariance 2 is not only relevant in 
the discussion of PDFs, but also for k± -factorization in 
general as well as for the cross sections which are k±- 
factorized [28-31] into an off-shell (fc^-dependent) par- 
tonic cross section it(| , k\ ) and a fc^-unintegrated parton 
density function 3 T{z 1 k 2 _): 

a = J ^-d 2 k ± a(^,k 2 ± )f{z,ki) (8) 

Here the partons generating a QCD hard process are off- 
mass shell. On-shell amplitudes in, say, dimensional rcg- 
ularization are supposed to be gauge- invariant, if not yet 
physical. The ensuing factorization of mass singularities, 
introduces a scheme- and perhaps gauge-dependence, to 
be canceled by (integrated) PDFs in a physical process in- 
volving hadrons. Thus, any gauge-dependence introduced 
in the PDFs is in a sense an artifact of the factorization 
procedure. 

A single off-shell gluon is not gauge-invariant. How- 
ever, experience with string theory [32] suggests that high- 
energy factorization could be a way of defining a phys- 
ical off-shell continuation, as residue of the Regge pole 
exchanged at the given (off-shell) momentum transfer. In 
such a case the k± cannot be assigned to a single gluon 
(except in some approximation) because the Reggeon is a 
compound state of (infinitely) many partons. Therefore, 
implementing such an idea in a formal definition is hard 
and further complicated by the fact that gluon Reggeiza- 
tion is infrared singular. 

The work on fcj_-factorization by [30] provided a gauge- 
invariant definition of off-shell matrix elements, based on 

2 this section is based on some remarks by M. Ciafaloni (un- 
published) during a discussion with John Collins and Yuri Dok- 
shitzer 

3 We use the classification scheme introduced in Ref.[l]: 
xQ(x,k\) describes DGLAP type unintegrated gluon distri- 
butions, xT{x,k 2 L ) is used for pure BFKL and xA(x, k\ , if) 
stands for a CCFM type or any other type having two scales 
involved. 



the Regge-gluon factorization idea. 4 The gluon Green 
function (and related uPDF) was defined so as to satisfy 
the (gauge-independent) BFKL equation, and the empha- 
sis was on defining the corresponding off-shell matrix ele- 
ments, given the physical cross-section. 

At leading-log level, CCH [30] noticed that the LO 
off-shell matrix elements could be defined by the high- 
energy limit of an on-shell six-point function (or an eight- 
point function in the two- k± case) whose expression was 
worked out in a physical gauge first, and then translated 
to the Feynman gauge. Because of their definition, the 
LO matrix elements are gauge-invariant and positive def- 
inite. At next-to- leading parton level Ciafaloni [33], and 
Ciafaloni and Colferai [34] noticed that one has to sub- 
tract, however, the leading kernel contribution (including 
gluon Reggeization) in order to avoid a rapidity diver- 
gence related to the ln(s) term in the total cross-section. 
This subtraction introduces a factorization-scheme depen- 
dence, mostly on the choice of the scale of the process, but 
not a gauge dependence. Recently the DESY group [35- 
37] went a long way towards completing this approach 
for DIS and jet production. The NLO matrix elements so 
defined are gauge- invariant, while positivity has not yet 
been thoroughly investigated, and is not guaranteed, be- 
cause of the subtraction. The latter is devised so as to put 
the whole energy dependence in the Green function. 

The CCFM equation employs a definition of uninte- 
grated density as a sum over physical final states, re- 
stricted to some angular region via angular ordering. This 
definition with rapidity cutoff is consistent with the sub- 
sequent analysis of matrix elements because the latter 
roughly provide upper and lower bounds on the rapid- 
ity integration, due to the angular coherence property. 
However, the relation of CCFM to BFKL was worked out 
at leading log (LL) only, and no complete attempt has 
been made so far to match this definition to exact next- 
to-leading log (NLL) calculations. In this case the energy 
dependence of the physical cross-section is shared between 
density and matrix elements, depending on the choice of 
the cutoff. 

The conclusion of the above considerations is, that 
any prediction for a physical process must be, obviously, 
gauge-invariant, however (unintegrated) PDF's are not 
guaranteed to be so. The formulation of fc^-factorization 
was meant to be gauge- invariant, and has been carried 
through at LL by [30,31] and at NLL level by [34,38,39]. 
It is not yet clear, whether gauge invariance is restored 
beyond that level. However, gauge-dependent definitions 
of PDF's with the corresponding matrix elements can be 
conceived also, provided their convolution reproduces the 
same (physical) cross-section. 



4 Note, that there are certain issues on gauge invariance 
which the authors of [30] and [24] have not been able to re- 
solve completely, but which will be a topic of future work. 
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3 Theoretical applications 

3.1 Comparison of available parameterizations 

The original CCFM splitting function given in eq.(l) in- 
cludes only the singular terms as well as a simplified treat- 
ment of the scale in a s , i.e. fcj_was used as the scale in 
the l/z term and the non-Sudakov factor whereas pj_was 
used in the 1/(1 — z) term and in the Sudakov form factor. 
Due to the angular ordering a kind of random walk in the 
propagator gluon fcj_can be performed, and therefore care 
has to be taken for small values of k±. Even during the 
evolution the non-perturbative region can be entered for 
k± < fcj"*. In the region of small k±, a s and the parton 
density are large, and collective phenomena, like gluon re- 
combination or saturation might play a role. Thus, the 
fast increase of the parton density and the cross section 
is tamed. However, for the calculation of the unintegrated 
gluon density presented here, a simplified but practical 
approach is taken: no emissions are allowed for k± < fc™* 
and q± < Qo- The limitation of fcjjs necessary for the 
calculation of the non-Sudakov form factor A ns in Eq.(3) 
and it ensures a finite value of a s (k±). Different choices 
of fc™* are discussed below. 

Following the arguments in [1], the scale in a s was 
changed to p± = q(l — z) everywhere, and the CCFM 
splitting function was extended to include also the non- 
singular terms [40,41]. The unintegrated gluon density 
at any x, k± and scale q is obtained by evolving nu- 
merically [6] a starting gluon distribution from the scale 
Qo according to CCFM to the scale q. The normaliza- 
tion ./V of the input distribution as well as the starting 
scale Qo, which also acts as a collinear cutoff to define 
Zmax = 1 — Qo/<l, need to be specified. These parameters 
were fitted such that the structure function F 2 as mea- 
sured at HI [42,43] and ZEUS [44,45] can be described 
after convolution with the off-shell matrix element in the 
region of x < 5 • 1CT 3 and Q 2 > 4.5 GeV 2 . Using 248 
data points a x 2 /ndf = 4.8, 1.29, 1.18, 1.83 for JS, J2003 
set 1,2,3 , respectively, is obtained. The following sets of 
CCFM unintegrated gluon densities are obtained: 

• JS (Jung, Salam [6]) 

The splitting function P gg of eq.(l) is used. The soft 
region is defined by fe™* = 0.25 GeV. 

• J2003 set 1 (Jung [41]) 

The splitting function P gg of eq.(l) is used, with fc™* — 
Qo fitted to kf 1 = Qo = 1.33 GeV. 

• J2003 set 2 

The CCFM splitting function containing also the non 
singular terms is used: 



10 GeV 



P(z,q, k)=a s (k\) ■ 
'(1-z) 



(9) 



+ z(l-z)/2)A ns {z,q,k) 



((i - *)V 



+ z{l-z)/2 



The Sudakov and non-Sudakov form factors were 
changed accordingly. The collinear cut is fitted to 
Qo = kf l = 1.18 GeV. 
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Fig. 2. Comparison of the different sets of unintegrated gluon 
densities obtained from the CCFM evolution as described in 
the text. In (a — c) the unintegrated gluon density is shown as 
a function of x for different values of k±at a scale of q = 10 



GeV. In (d - f) the ratio R = x ^ (x ^ ±, l ) 
x for different values of k±is shown. 



as a function of 



• J2003 set 3 

CCFM splitting function containing only singular 
terms but the scale in a s is changed from k± to 
p± for the l/z term. The collinear cut is fitted to 
Qo = fcj"* = 1-35 GeV. The problematic region in 
the non-Sudakov form factor in cq.(3) is avoided by 
fixing a s {[i r ) for \i r < 0.9 GeV. 

A comparison of the different sets of CCFM unintegrated 
gluon densities is shown in Fig. 2. It is clearly seen, that 
the treatment of the soft region, defined by k± < fe™* 
influences the behavior at small x and small k±. 

Also the LDC model describes F 2 satisfactorily, but the 
corresponding unintegrated gluon densities are somewhat 
different. One major difference as compared to CCFM is 
that LDC can also include quarks in the evolution, and can 
therefore also reproduce Fi in the valence region of high 
x. In Fig. 3 three different unintegrated gluon densities for 
the LDC approach are presented. The standard set refers 
to the full LDC including quarks in the evolution and the 
full gluon splitting function, whereas for the gluonic set 
and the leading set only gluon evolution is considered with 
only singular terms in the splitting function for the latter. 
All three alternatives have been individually fitted to F 2 
in the region x < 0.3, Q 2 > 1.5 GeV 2 for standard and x < 
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Fig. 3. Comparison of the different sets of unintegrated gluon 
densities obtained within LDC at scale of q = 10 GeV. Stan- 
dard refers to the full LDC including quarks in the evolution 
and the full gluon splitting function. For gluonic and leading 
only gluon evolution is considered with only singular terms in 
the splitting function for the latter. Also shown is the J2003 
set 1 for comparison (divided by ir). 



0. 013 and Q 2 > 3.5 GeV 2 for gluonic and leading. In LDC 
there is only one relevant infrared cutoff, fc^o, which limits 
the p± of emitted gluons. This has been fitted to 0.99, 1.80 
and 1.95 for standard, gluonic and leading respectively. 
No cut on the transverse momenta of the virtual gluons 
is applied and the argument /Lt r in a s is set to /i r = p± 
which is then always larger than the cutoff k± • 

3.2 Semi analytical insight into the CCFM equation 

The CCFM equation interlocks in a rather complicated 
way the two relevant scales, i.e. the transverse momentum 
k± of the parton and the hard scale q which is related 
to the maximal emission angle. Due to this complexity 
the existing analyses of the CCFM equation are based 
upon numerical solutions. After performing some approx- 
imations it is however possible to obtain semi analytical 
insight into the CCFM equation and we would like to con- 
sider the following two cases: 

1. The single loop approximation (SLA) [46-48], which 
corresponds to the DGLAP limit. 

2. The CCFM equation at small x with consistency con- 
straint (CC) [3,49]. 

3.2.1 The single loop approximation 

This approximation corresponds to setting the non- 
Sudakov form-factor equal to unity and to the replace- 
ment of the angular ordering by a q± ordering. It can 
be a reasonable approximation for large and moderately 
small values of x. In Fig. 4 we show a comparison of the 



unintegrated gluon density obtained in the full CCFM 
and the single-loop approximation (using the same input 
distributions). 

In SLA the CCFM equation can be partially diagonalized 
by the Fourier-Bessel transform of the unintegrated gluon 
distribution A(x, k±,q) [50] 

POO 

A(x,kj_,q) = / dbbJ f) {k^,b)A{x,b,q) (10) 
Jo 

POO 

A(x,b,q)= / dkj_kj_J a (kj_,b)A(x,k±,q) (11) 
Jo 

with b being the impact parameter and the integrated 
gluon distribution is given by: 

xg(x,Q 2 ) = 2A(x,b = 0,q = VO 1 ) (12) 

The transverse coordinate representation also partially di- 
agonalizes the CCFM equation extended for polarized un- 
integrated parton distributions as will be discussed in a 
separate section below. Due to the absence of the 1/z 
term in the polarized splitting function it will be possi- 
ble to utilize this representation beyond the single loop 
case (see section 2.6). The tranverse coordinate represen- 
tation has also been used for the analysis of the CCFM 
equation in SLA for the unintegrated gluon distributions 
in a photon [51]. 

The CCFM equation in SLA takes the following form 
in the b representation: 

A(x, b ,q) = A° { x,b) + f 2 ^^ f^zPM- 

Jql Q 27T J Q Z 

[0{z-x)J Q [b{l- z)q]A( K - z ,b,q) -zA(x,b,q)} (13) 

where for simplicity we have neglected the quark contri- 
bution. At b = this equation reduces to the conventional 
DGLAP evolution equation in the integral form. Equa- 
tion (13) can be solved in a closed form using the moment 
function 

f u (b,q) = ( dxx^AixAq) (14) 
Jo 
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Different approximations of eq.(13) are related to for- 
malisms used e.g. in studies of Drell-Yan pairs, and can 
give more insight into the properties of the solutions to 
the CCFM equation. Thus approximations in the Bessel 
function and the Sudakov form-factor gives the relation 
[50]: 



a, , T g (b=l/k ± ,q) n-^/t 

A{x,k ± ,q) ~ -2- 



Jo 



dzPgg(z)- 



2tt ' z \z 

where the Sudakov form- factor T g (b,q) is defined by: 



T g (b,q) = exp 



« 2 dq^astf) ^-VCM 
1/6 2 q 2 2w 



dZZPgg{z) 



(16) 

Neglecting also contributions to the Sudakov form-factor 
for large q 2 gives: 



A(x, k±,q) ~ 



d[T g {b=l/k ± ,g)xg{x,kl)] 
dk\ 



(17) 



The approximate expressions (17) and (15) are similar to 
those discussed in [25,52,23]. It turns out that expression 
(15) gives a reasonable approximation of the exact solution 
of the CCFM equation in SLA while expression (17) can 
generate negative result for large k± and large x ~ 0.1 
[50]. 



3.2.2 CCFM equation with consistency constraint 

We shall consider now the CCFM equation in the small 
x limit keeping only the singular l/z part of the split- 
ting function P gg (z) and neglecting the Sudakov form- 
factor. We shall also impose the consistency constraint [3, 
49] which is known to generate the dominant part of the 
sub-leading BFKL corrections. The integration limit (s) in 
the CCFM equation are now constrained by the following 
competing conditions: 5 

1. Angular ordering (AO) <-> Zi-\q~i-\ < q~i. 

2. Consistency constraint (CC) ^> q 2 < kj_ t /zi. 

It can easily be observed that CC takes over AO for k 2 ^ < 
f/z. 

The structure of the CCFM equation at small x with 
CC is different in the regions k± < q and k± > q. At k± < 
q the unintegrated distribution A(x, k±,q) is independent 
of q, i.e.: 

A(x,k±,q) -+T(x,k±) (18) 

after adopting the leading ln 2 (fc^/g 2 ) approximation of 
the Sudakov form factor, while for k± > q we get the 



5 The single loop approximation is extended, since the q or- 
dering is replaced by angular ordering 



following expression after adopting the leading double 
In 2 {k\/q 2 ) approximation: 6 



A(x, kj_, q) = T(x, kj_) exp 



3q s 
" 2tt 



\n\ki/q 2 ) 



(19) 



where for simplicity we keep fixed a s . The single scale 
function !F(x,k±) satisfies the BFKL-likc equation with 
sub-leading corrections. We found in this way that impos- 
ing the consistency constraint and the double ln 2 (fc 2 _/g 2 ) 
approximation in the region k±_ > q we reduce the two- 
scale problem to the single-scale one and to the BFKL-like 
dynamics. The novel feature of the CCFM framework is 
however the exponential suppression of the unintegrated 
distribution in the region k± 3> q (cf. equation (19)) due 
to the double ln(kl/q 2 ) effects [15]. They are of course 
formally sub-leading at small x. 



3.3 Effects of phase space constraints in BFKL 

The leading logarithmic (LL) BFKL formalism resums 
terms in the perturbative series of the form (a s ln(s/so)) n , 
where s is the square of the center of mass energy for the 
hard scattering and sq some perturbative scale separating 
the evolution of the i-channel exchange from the matrix 
elements. These logarithms arise due to the emission of 
gluons from the i-channel exchange. For the scattering 
of two particles paPb — » k a kbki where fcj are the mo- 
menta of the gluons emitted from the BFKL evolution, 
we have s = 2paPb and sq is often chosen as sq = k a j_kb± 
with k a ± (kb±) the transverse part of k a {kb respectively). 
In deep inelastic scattering large s corresponds to small 
x of the probed parton. For hadronic dijet production, 
large s corresponds to large separation in rapidity between 
the leading jets, and therefore to moderate values of x, 
where normal DGLAP evolution of the partons is valid 
(and therefore the standard PDFs can be used [53]). The 
next -to-leading logarithmic corrections [54,33] consist of 
terms proportional to a s (a s ln(s/so))™, i.e. suppressed by 
one power of a s compared to the LL component. The log- 
arithms resummed in the BFKL approach correspond to 
the enhanced terms in scattering processes for large center 
of mass energies and also the enhanced terms in the de- 
scription of the small x behavior of the gluon distribution 
function. 

When confronting BFKL predictions with data, sev- 
eral points are worth observing. First of all, present day 
colliders do not operate at "asymptotic energies" where 
the high energy exponent dominates the BFKL prediction 
under the assumption that the coupling can be held fixed 
and small, leading to a prediction of an exponential rise in 
cross section with an intercept of a s 41n2, with a s = ^Si. 
For example at HERA, the separation between the struck 
quark and the forward jet can reach up to about four units 
of rapidity, whereas the measurable jet separation at the 



6 Double ln(fc|/(j 2 ) terms appear in the Sudakov form fac- 
tor for exclusive cross sections. They are not present in the 
inclusive cross section, which is oc F2 
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Tevatron is up to six units. This is not asymptotically 
large. Secondly, it should be remembered that the loga- 
rithms rcsummed are kinematically generated, and in the 
derivation of the standard analytic solution to the BFKL 
equation, the transverse momentum of the gluons emitted 
from the BFKL evolution has been integrated to infin- 
ity. It is therefore apparent that any limits on the phase 
space probed in an experiment can have a crucial impact 
on the theoretical prediction. Such limits can either be 
the cuts implemented in the measurement or the limits 
on the available energy at a collider. The total available 
energy will affect the impact factors, while taking into 
account also detailed energy-momentum conservation in 
each gluon emission will in addition affect the BFKL ex- 
ponent. Taking hadronic dijet production as an example, 
the energy constraint will obviously not just limit the pos- 
sible rapidity separation of the leading dijets, but also the 
amount of possible radiation from the BFKL evolution, es- 
pecially when the leading dijets are close to the kinemati- 
cal boundary. For a multi-particle final state described by 
two leading dijets with transverse momentum and rapid- 
ity {p a /b±,y a /b) and n gluons described by (fcu,J/»), the 
total energy of the event is given by s = x a XhS where s is 
the square of the total energy of the hadron collider and 



Pa± Va , \ ^ «i_L Vi . Pb± 



(20) 



While it can be argued that the contribution to s from 
the gluons emitted from the BFKL evolution is sublead- 
ing compared to the contribution from the leading di- 
jets, it is not obvious that the effect on the cross section 
is small, simply for the reasons mentioned above: ignor- 
ing the contribution to the parton momentum fractions 
will resum logarithmically enhanced contributions from 
regions of phase space that lie outside what can be probed 
at a given collider. Here we have taken as an example dijet 
production at a hadron collider, but a similar effect will be 
found for any BFKL evolution, whether it describes 7*7*, 
ep, or pp physics. 

The iterative approach of Ref. [55,56] to solving the 
BFKL equation at leading logarithmic accuracy allows not 
only a study on the partonic level of BFKL in processes 
with complicated impact factors, where it might be diffi- 
cult to get analytic expressions for the cross section. The 
method also allows for the reconstruction of the full final 
state configurations contributing to the BFKL evolution, 
and therefore it is possible to study quantities such as mul- 
tiplicities and distribution in transverse momentum of the 
emitted gluons [57] . Only this reconstruction of the full fi- 
nal state allows for the observation of energy and momen- 
tum conservation. The effects of energy and momentum 
conservation have been studied in several processes [58- 
60] . When no phase space constraints are imposed, the it- 
erative solution reproduces the known analytic solution to 
the BFKL equation. This iterative approach has recently 
been generalized [61,62] to solving the NLL BFKL equa- 



tion thereby joining other approaches [63-67] in studying 
effects of NLL corrections. 

The effects on the total center of mass energy of consid- 
ering the full multi-gluon BFKL final state in gluon-gluon 
scattering is seen in Fig. 5. We have plotted the result of 
considering only the two leading dijets (i.e. ignoring the 
sum in Eq. (20)), and from considering the full BFKL fi- 
nal state (i.e. using the full expression in Eq. (20)) (see 
Ref. [57] for more details). Fig. 5 shows the average en- 
ergy for a BFKL dijet event as a function of the rapidity 
separation of the leading dijets, when the BFKL gluon 
phase space is unconstrained. The standard analytic solu- 
tion to the BFKL equation implicitly assumes that all of 
this phase space is available. It is clear from Fig. 5 that 
the energy taken up by BFKL radiation is significant com- 
pared to the center of mass energy at present and planned 
colliders. For example at four units of rapidity, which is 
the upper limit of HERA at present, the average energy 
of a BFKL dijet event is about y/§ = 1 TeV according to 
Fig. 5, which is far beyond the maximum energy available. 
At the HERA center-of-mass energy \/§ = 300 GeV the 
rapidity range would be less than one unit and thus leave 
very little phase space for additional emissions. Therefore, 
any constraint on the BFKL radiation from e.g. overall 
energy conservation will have an impact on BFKL phe- 
nomenology predictions at such colliders. In fact, it is 
found that if energy and momentum conservation is satis- 
fied, by using the full Eq. (20) when calculating hadronic 
dijet production at the LHC, then the exponential rise 
in cross section as a function of the rapidity separation 
found for gluon-gluon scattering (when the BFKL gluon 
phase space is integrated to infinity) is moderated to an 
almost no-change situation compared to the fixed lead- 
ing order QCD prediction. Other BFKL signatures, like 




Fig. 5. The average center of mass energy in gg — > gg scatter- 
ing with (red/ dashed) and without (green/ dotted) BFKL evo- 
lution of the t channel gluon, with px m i n = 20 GeV for the 
dijets and a s — 0.1635. Also plotted is the hadronic center 
of mass energy squared for the Tevatron ((l.8TeV) 2 ) and the 
LHC((14TeV) 2 ). 
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3.4 The saturation scale 
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Fig. 6. The 0-, 1, 2-, 3- and 4-jet parton-level cross sections 
as a function of the cutoff hr, for a rapidity span of Ay = 5 
and p_L m i n = 20 GeV for the leading dijets. Also shown is the 
analytic 0-jet prediction valid for small hr. 



the increasing dijet angular de-correlation with increasing 
rapidity separation [68-70], are still present. 

With the iterative method of solving the LL BFKL 
equation it is of course also possible to calculate jet rates 
and transverse momentum distributions (since full infor- 
mation of the final state configuration is obtained) arising 
from the BFKL dynamics. Below we present results on 
the jet rates in gluonic dijet production (on the partonic 
level, i.e. with the full gluonic phase space assumed in the 
standard analytic solution of the BFKL equation) with 
a BFKL chain spanning 5 units of rapidity using a very 
simple jet definition. We simply let any gluon with a trans- 
verse momentum greater than some cut off fin define a jet 
(this is a reasonable jet definition since at leading loga- 
rithmic accuracy the emitted gluons are well separated in 
rapidity). These jet rate predictions will change, once the 
partonic cross section is convoluted with parton density 
functions. The jet rates of Fig. 6 are the ones responsible 
for the increase in the center of mass energy of a BFKL 
event over a simple LO configuration seen in Fig. 5, but 
they are also responsible for the rise in cross section pre- 
dicted from the BFKL dynamics when the BFKL gluonic 
emission is unbounded. In Ref. [57] it is found that for 
gluon-gluon scattering (with a minimum transverse mo- 
mentum of the leading dijets of 20 GeV) with a BFKL 
exchange, one can expect a BFKL gluon emission density 
of about one hard (k±i > 20 GeV) gluon for every two 
units of rapidity spanned by the BFKL evolution, when 
the energy of the event is unconstrained. This amount of 
radiation is implicitly assumed in the standard analytic 
solution of the LL BFKL equation. 

In conclusion, carefully taking energy-momentum con- 
servation into account dramatically modifies the strong 
increase for small x, predicted in the leading log BFKL 
approach. 



The parton saturation idea is realized with the help of a 
nonlinear evolution equation in which the gluon splitting 
is described by a linear term while the negative nonlinear 
term results from the competing gluon recombination (see 
also [28]). The Balitsky-Kovchegov (BK) equation [71,72] 
was derived for deep inelastic scattering of a virtual pho- 
ton on a large nucleus by the resummation of multiple 
pomeron exchanges in the leading logarithmic approxi- 
mation (when a s ln(l/:r) ~ 1 and x ~ Q 2 /s) in the large 
iV c limit. It is an equation for the dipole-proton forward 
scattering amplitude iV(x, y, Y) where x, y are the end 
points of the qq dipole and Y = In 1/x is the rapidity of 
the process. The BK equation has the following integro- 
differential form [72] 



cW(x,y, Y) 

dY 
N(y,z,Y) 



d 2 z(x 



[iV(x,z,V)- 



(x-z) 2 (y-z) 2 
- iV(x, y, Y) - iV(x, z, Y)N{y, z, Y)} (21) 



where a s = 3a s /ir and is fixed in the leading lnl/a; ap- 
proximation. The linear term in (21) is the dipole version 
of the BFKL equation whereas the quadratic term de- 
scribes the gluon recombination. Instead of x, y one often 
uses their linear combinations: r = x — y which is the 
size of the dipole and b = h (x + y) the impact param- 
eter. Equation (21) can be easily solved when using the 
approximation of the infinitely large nucleus, i.e. assum- 
ing that the amplitude N(x,y, Y) = N{\r\ : Y) depends 
only on the size of the dipole but not on the impact pa- 
rameter [73-75] . The more complicated case with the full 
impact parameter dependence has been analyzed recently 
[76]. For the 6-independent and cylindrically symmetric 
solution, N(r,Y) — N(r,Y), Eq. (21) can be rewritten in 
momentum space in a much simpler form after performing 
the following Fourier transform 



(j){k,Y) = j — exp(-ik-r 



N(r, Y) 



dr 



J (kr)N(r,Y) 



(22) 



where Jo is the Bessel function. In this case the following 
equation is obtained 



d</>(k,Y) 
dY 



a s (K®4>)(k,Y) 



</> 2 (fc,n (23) 



Here expression (K £g) <fi)(k,Y) means the action of the 
usual BFKL kernel in the momentum space onto the func- 
tion (j>(k,Y). Let us briefly analyze the basic features of 
N(r,Y) and 4>(k, Y). In Fig. 7 we plot the amplitude 
N(r, Y) as a function of the dipole size r for different 
values of rapidity Y. The amplitude N(r,Y) is small for 
small values of the dipole size. It is governed in this regime 
by the linear part of equation (21). For larger values 
of dipole sizes, r > 1/Q S (Y) the amplitude grows and 
saturates eventually to 1. This is the regime where the 
nonlinear effects are important. As it is clear from the 
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Fig. 7 the saturation scale grows with rapidity Q S (Y). It 
means that with increasing rapidities the saturation oc- 
curs for smaller dipoles. It has been shown that the growth 
of the saturation momentum is exponential in rapidity 
Q s {y) = Qo exp(AF) with A ~ 2a s being a universal co- 
efficient and governed by the equation. The normalization 
Qo on the other hand is dependent on the initial condition 
N(r, Y = 0). We note however, that when the rapidity Y 
is not too large the initial conditions are still important. 
In this region the coefficient A can still depend on the 
rapidity [74]. 

It has to be stressed that in the leading-log x approxi- 
mation the strong coupling constant is fixed. The running 
of the coupling, although being a next-to-leading-log x ef- 
fect, is obviously more physical. In this case the rapidity 
dependence of the saturation scale is changed. We adopt a 
natural approximation that the local exponent of the sat- 
uration scale A(Y) = d\n(Q s (Y) / A) jdY takes the form 
X(Y) = 2a s (Q 2 s (Y)) where A = A QC d- The above form is 
motivated by the leading logarithmic result with the fixed 
coupling as discussed before, i.e. Q S (Y) = Qoexp(AF). 

Thus, we have 



d\n(Q s {Y)/A) 
dY 



12 



boHQ s {Y)/A) 



(24) 



with the initial condition Q s (Yq) — Qo and Yq chosen in 
the region where scaling sets in. The solution takes the 
form 



•I "M>( \/|V- Y o) 



Y > Yn 



(25) 

where Lo = ln(Qo/^)- Thus, the exponential dependence 
on the rapidity Y of the saturation scale is changed in 
the running coupling case to the milder behavior [75]. In 
the phenomenological analysis of the HERA data based 
on the BK equation [77], the running of the a s has been 
included. It has been shown that, in the limited range 
of x, the effective power governing the behavior of the 
saturation scale can be still fitted using the simple form 
Qs(Y) ~ exp(AoY) with Ao — 0.18, which is close to the 
value used in the Golec-Biernat and Wiisthoff saturation 
model [78,79]. 

In what follows we consider the solution to the 
Balitsky-Kovchegov equation in the leading-log x fixed 
coupling case. In Fig. 8 we plot the solution k(f>(k,Y) in 
momentum space as a function of k for increasing rapidi- 
ties. We have imposed the initial condition of the form of 
the delta function, localized at some k = k Q = 1 GeV. 
We compare the solution of the BK equation with that 
of the linear BFKL which exhibits the strong unlimited 
diffusion into both infrared and ultraviolet regimes. The 
BFKL solution is symmetric and peaked around the initial 
value k = fco- The solution to the BK equation shows on 
the other hand a suppression of the diffusion into the low 
momenta. We clearly see that the peak of the distribution 
moves with increasing rapidity to the higher momenta. 
One can identify the value of the momentum k at which 
the maximum occurs as the saturation scale Q S (Y). At 



large momenta k where the nonlinearity in the BK equa- 
tion does not play a role, the two solutions BFKL and BK 
are close to each other. The overall height of the distribu- 
tion is strongly damped with respect to the linear case. 




Fig. 7. The amplitude N(r, Y) as a function of the dipole size 
r for two different rapidities Y = 1,5 



The solution to the BK equation shows also another inter- 
esting feature, namely the property of scaling (e.g. [80,75, 
81]). In the nonlinear regime r > 1/Q S (Y) the amplitude 
depends only on one combined variable instead of r and 
Y separately. 



N(r,Y) ee N(rQ a (Y)) 



(26) 



This is also property of the GBW saturation model, 
though in the latter case the scaling was present for all 
values of r. In the case of BK equation, scaling only oc- 
curs in the saturation domain that is for large values of 
dipole sizes. 



3.5 Non-linear evolution versus HERA data 

A new approach to a global QCD analysis based on the 
non-linear QCD evolution by Balitsky-Kovchegov (BK) is 
presented in [74,77]. The BK equation improved by the 
DGLAP corrections for small dipole sizes (independently 
of impact parameter) is in fact very successful in describ- 
ing the low x part of the structure function F2 at HERA. 
In the following a brief summary of the results in Ref. [77] 
is given. 

With the initial conditions specified at xo = 10~ 2 the 
BK equation (without impact parameter dependence) is 
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Fig. 9. Fit to the F2 structure function. The dashed line is a result obtained without the DGLAP corrections. The data are from 
Refs. [82,45,83,84] 



solved numerically towards smaller x. The impact param- 
eter dependence is restored using a rescattering ansatz of 
the Glauber-type. All existing low x (x < 0.01) data on 
the F2 structure function are reproduced with resulting 
X 2 /ndf ~ 1 (Fig. 9). Only two parameters and a few 
fixed ones (associated with the initial conditions) are used 
for the fit. The fitted parameters are the effective proton 
radius, entering the Gaussian impact parameter distribu- 
tion, and the scale at which the DGLAP corrections are 
switched on (0(1 GeV)). The DGLAP corrections are im- 
portant for large photon virtualities only and reach up to 
15%. The low Q 2 (of the order of a few GeV 2 and below) 
data are described solely by the BK equation. 



The main fitting parameter used for the fit is an ef- 
fective proton radius R. The optimal fit is achieved at 
R ~ 0.3 fm, the radius which is much smaller than the 
electro-magnetic radius of proton. On one hand, this small 
proton radius might be an artifact of the approximations 
used. On the other hand, it may indeed indicate a small 
size dense gluon spot inside proton. Such scenario arises in 
several other models for high energy scattering off proton 
[85,86]. 

The approach based on non-linear QCD evolution al- 
lows the extrapolation of the parton distributions to very 
high energies available at the LHC as well as very low 
photon virtualities, Q 2 < 1 GeV 2 . 



In DIS the Pomeron intercept is obtained by a mea- 
surement of A = din F2 / dhi(l / x) . For large photon vir- 
tualities the fit based on the BK equation reproduces the 
HERA data with A ~ 0.3 - 0.4, the hard BFKL inter- 
cept. In the small Q 2 region the non-linear terms in the BK 
equation are reflected in the smaller A values at smaller 
Q 2 . Fig. 10 presents a prediction for A at at smaller Q 2 and 

smaller x. Fig. 11 presents results for x ~ 10~ 4 . In this re- 
gion A decreases strongly for small Q 2 , but varies relatively 
slowly with x. In fig. 10 we see, however, that for smaller 
values of a;, A decreases more strongly with x, for fixed Q 2 , 
tending to zero in agreement with the unitarity constrain. 
At Q 2 well below lGeV 2 and x ~ 10~ 6 , A ~ 0.08 - 0.1. 
This value of A coincides with the "soft Pomeron" inter- 
cept. Thus the nonlinear evolution provides a solution to 
the problem of hard-soft Pomeron transition. 



3.6 Multiple Interactions in non-ordered cascades 

At high energies the perturbative jet cross section in pp 
collisions becomes larger than the total cross section. This 
implies that there are often several hard sub-collisions in 
a single event. Therefore correlations become important, 
and the observed "pedestal effect" implies that the hard 
sub-collisions are not independent [87], indicating an im- 
pact parameter dependence such that central collisions 
have many mini-jets, while peripheral collisions have fewer 
mini-jets [87]. Also at HERA the final state properties in 
photoproduction cannot be reproduced without assuming 
multiple hard scattering [88-90] . At higher Q 2 the indica- 
tions for multiple scattering are reduced, and thus HERA 
offers a unique possibility to study how, with decreasing 
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Fig. 8. The Fourier transform of the amplitude N(r, Y) as a 
function of the momentum k for different values of rapidities Y 
increasing from 1 to 10. Dashed lines correspond to the solution 
to the BFKL equation whereas solid lines to the full nonlinear 
BK equation. 



Q 2 , multiple interactions become more and more impor- 
tant, until eventually a situation similar to pp collisions is 
reached for Q 2 = 0. 

In a non-fcj_ -ordered BFKL ladder, it is possible to 
have two (or more) local fc^-maxima, which then corre- 
spond to two different hard sub-collisions. Thus there are 
two different sources for multiple interactions: It is possi- 
ble to have two hard scatterings in the same chain, and 
there may be more than one chain in a single event. The 
BFKL, CCFM or LDC formalism can be used to estimate 
multiple collisions in a single chain. The symmetric prop- 
erties of the LDC model for DIS makes it especially suited 
to be applied to pp collisions, and in Ref. [91] it is demon- 
strated that it is possible to deduce the average number 
of chains in pp scattering from data on deep inelastic ep 
scattering. 

The LDC model can, however, not determine the corre- 
lations between the chains. Uncorrelated chains would be 
described by a Poissonian distribution, but the observed 
pedestal effect, mentioned above, makes it more likely 
that central collisions have more, and peripheral collisions 
fewer, chains. The analysis by Sjostrand and von Zijl [87] 
favors an impact parameter dependence described by a 
double Gaussian distribution. It turns out that this distri- 
bution leads to a geometric distribution in the number of 
sub-collisions, with the tail suppressed by energy conser- 
vation. Some predictions for mini-jet multiplicity and the 
pedestal effect in pp collisions are presented in Ref. [87], 
assuming such a geometric distribution for the number of 
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chains in a single pp event. Further work is in progress, and 
it would be very interesting to test these ideas, not only 
in pp or pp collisions, but also in ep scattering, varying Q 2 
from the DIS region to photoproduction. 
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3.7 Spin dependent unintegrated parton distributions 

The basic, universal quantities which describe the inclu- 
sive cross-sections of hard processes within the QCD im- 
proved parton model are the scale dependent parton dis- 
tributions. These parton distributions or distribution am- 
plitudes describe how the momentum of the nucleon is 
distributed among its constituents, i.e. quarks and gluons. 

Polarized parton distributions are a probabilistic mea- 
sure for the distribution of the nucleon's longitudinal spin 
(helicity) among its constituents. More precisely, one de- 
fines polarized parton distributions as the difference of the 
probability density to find a parton / with its longitudi- 
nal spin parallel aligned minus the probability density to 
find the same parton with its longitudinal spin antiparallel 
aligned relative to the spin of the nucleon: 

Af = /tt - fn • (27) 

These parton distributions conventionally only depend on 
x and Q 2 , but just as for the spin-independent case it may 
be beneficial to also consider fc^-unintegrated polarized 
parton distributions. An evolution equation analogous to 
CCFM has been derived for the unintegrated gluon dis- 
tribution in [92] and the result is quite similar, although 
contrary to the unpolarized case there is no non-Sudakov 
form factor since the polarized splitting function does not 
have al/z pole. 

One can show that the principles discussed in [92] with 
slight modifications also apply for the case of including 
quarks in the evolution. Thus one arrives at a complete 
set of evolution equations along the lines of CCFM for 
the unpolarized case. The CCFM evolution for polarized 
gluons is called pCCFM evolution equation [93]. 

It can easily be shown that in the small x limit the 
pCCFM equation formulated in [92] generate the double 7 
ln 2 (l/x) for distributions integrated over transverse mo- 
mentum of the partons. Their detailed structure is how- 
ever different from the collinear QCD expectations [94,95]. 
One can modify the pCCFM equations in order to incorpo- 
rate those expectations and make contact with the evolu- 
tion equations in the integrated case containing Altarclli- 
Parisi + ladder contributions which have been discussed in 
[94] . These modifications contain the following steps [93] : 

1. In order to get the expected double logarithmic limit 
of the integrated distributions it is sufficient to replace 
the angular ordering constraint 0{Q — z|qj_|) by the 
stronger constraint 0(Q 2 — zq\) in the corresponding 
evolution equations for integrated distributions. 

2. The argument of a s will be set equal to q\ instead of 

g 2 = gi(l r *) 2 . 

3. The non-singular parts of the splitting function(s) will 
be included in the definition of the Sudakov form- 
factor^). 

4. Following Ref. [94] we include the complete splitting 
functions P a b{z) and not only their singular parts at 
z = 1 and constant contributions at z = 0. 

7 The double ln 2 (l/:r) terms come from non- ladder 
bremsstrahlung terms. 



5. We represent the splitting functions AP ab {z) as: 
AP ab (z) = AP ab (0) + AP ab {z) where AP o6 (0) = 0. 
Following [94] we shall multiply AP ab (Q) and AP ab (z) 
by 0(Q 2 —zq 2 _) and 0{Q 2 — qj_) respectively in the inte- 
grands of the corresponding integral equations. Follow- 
ing the terminology of Ref. [94] we call the correspond- 
ing contributions to the evolution kernels the 'ladder' 
and 'Altarelli-Parisi' contributions respectively. 

6. We shall 'unfold' the eikonal form factors in order to 
treat real emission and virtual correction terms on 
equal footing. 

After those modifications one arrives at an evolution equa- 
tion for the unintegrated polarized parton distributions 
which includes the complete LO Altarelli-Parisi and the 
double ln 2 (l/x) effect generated by ladder diagrams in a 
consistent way, i.e. if one integrates the evolved uninte- 
grated parton distributions over the transverse momen- 
tum fcj_the result will be the same as if one had done an 
evolution with the integrated parton distributions using 
'Altarclli-Parisi+laddcr' evolution equation. It means that 
the corresponding diagram between evolution and trans- 
verse momentum integration commutes. 

One can utilize the fact that the pCCFM equation can 
be (partially) diagonalized by the Fourier-Bessel trans- 
form. It turns out that the difference between the inte- 
grated and the unintegrated evolution equation in Fourier- 
space is only a single factor Jo(&_i_g_i_(l — z )), where b± is 
the transverse impact parameter conjugate to the trans- 
verse momentum of the parton, q± the transverse evolu- 
tion scale, z the momentum fraction and Jo the Bessel 
function of order 0. The evolution equation for the inte- 
grated case is simply restored by setting b± = 0. 

There is a third contribution to the evolution of unin- 
tegrated polarized parton distributions which is not cov- 
ered by the 'Altarelli-Parisi + ladder' approximation of 
the modified pCCFM equation, these are the non-ladder 
bremsstrahlung contributions. A general method of imple- 
menting the non-ladder bremsstrahlung corrections into 
the double logarithmic resummation was proposed by 
Kirschner and Lipatov [96,97]. For unintegrated polar- 
ized parton distributions they have been implemented in 
Ref. [94]. In the unintegrated case one can simply add 
them by analogy to the Altarelli-Parisi and ladder con- 
tribution by inserting the factor Jo(b±q±(l — z)) because 
then again one obtains perfect commutativity of the di- 
agram between the integrated and unintegrated parton 
distributions. 

Putting all three contributions together (Altarelli- 
Parisi, ladder + non-ladder) one obtains a set of linear 
integral equations for unintegrated polarized quark and 
gluon distributions. In Fig. 12 we show the evolution of 
the k± dependence for the triplet contribution : 

Ah = ^(Au + Au- Ad- Ad) . (28) 

The input distributions at Qq = 0.26 GcV 2 arc taken 
from [98, 99] and are compared to the evolved distributions 
at Qq = 10.0 GeV 2 . The width of the initial transverse 
momentum dependence a has been chosen to be 1 GeV. 



14 



Jeppe Andersen et al.: Small x Phenomenology - Summary and status 



triplet contribution, ^dependence 




Fig. 12. Transverse momentum dependence for the triplet con- 
tribution Af 3 = | (Au + Au - Ad- Ad) using pCCFM evo- 
lution including non-ladder contributions. The thin lines show 
the input distributions for Qq = 0.26 GeV 2 for x — 0.1 (solid) 
and x = 0.01 (dashed), while the bold lines show the same dis- 



tribution evolved to 



10 GeV^ 



For the simulation the Altarclli-Parisi, ladder and non- 
ladder contributions all have been included. It is seen that 
due to the evolution the kj_ dependence is broadening away 
from a Gaussian behavior to a more exponential decay. 



3.8 J/^-production and polarization effects 

In the following we consider J/ip meson production in ep 
deep inelastic scattering in the color singlet model using 
/c_i_-factorization. It should be noted that heavy quark and 
quarkonium cross section calculations within the collinear 
factorization in fixed order pQCD show a large discrep- 
ancy (by more than an order of magnitude) [100-103] to 
measurements at the Tevatron for hadroproduction of J / ip 
and T mesons. This fact has resulted in intensive theoret- 
ical investigations of such processes. In particular, it was 
suggested to add an additional transition mechanism from 
cc-pairs to J/ip mesons, the so-called color octet (CO) 
model [104], where a cc-pair is produced in a color octet 
state and transforms into the final color singlet (CS) state 
by the help of very soft gluon radiation. The CO model 
is based on the general principle of the non-relativistic 
QCD factorization (see [104]). The physics behind this 
factorization is simple: the heavy quark-antiquark meson 
is produced at distances which are not so short as the dis- 
tances for heavy quark-antiquark production (which are 
of the of l/2mg, where toq is the mass of heavy quark). 
Indeed, we can easily estimate that the typical distances 
for e.g. J/ip production is about l/a s (mg)mQ. These dis- 
tances are much longer than the distances of the typi- 
cal hard process but they are still much shorter than the 
hadronization distances. Therefore J/ip production is still 
under control of perturbative QCD but on the other hand 
it could be accompanied by a highly non-perturbative pro- 
duction of soft gluons. By adding the contribution from 
the CO model and fitting the free parameters one was able 



to describe the existing data on the production of J/ip pro- 
duction at the Tevatron. However, in recent years, we have 
seen a lot of difficulties of the CO model. The first and the 
most disturbing is the fact that the fit with the CO model 
gives values of wave functions at the origin which are in 
contradiction with the non-relativistic (NR) QCD hierar- 
chy where each non-perturbative CO matrix element has 
a definite order of magnitude in the relative heavy quark 
velocity. The qualitative prediction for the CO model is 
the transverse polarization of the produced J/ip since the 
main contribution of the CO model to J/ip production 
in pp-collisions comes from gluon-gluon fusion with trans- 
verse polarized gluons. The second important question is 
about the NR QCD factorization itself. Is the heavy quark 
mass really large enough to have well separated scales, 
l/2m<5 and l/2a s (mQ)mQ, or is a special selection needed 
as suggested in Rcf.[105] The CO model has been applied 
earlier in the analysis of inelastic J/ip production [106, 107] 
at HERA. However, as noted in Ref. [107], the results from 
Rcf. [106] and [107] do not agree with each other. Also the 
results obtained within the usual collinear approach and 
the CS model [108-111] underestimate the experimental 
data by factor about two. 

First attempts to investigate the J/ip polarization 
problem in ep- interactions at HERA and in pp- interactions 
at the Tevatron were made in [112-116] using the k±- 
factorization approach. An extensive analysis of the pro- 
duction of J/ip, Xc and T mesons (including the polar- 
ization properties) in pp-collisions has been recently pre- 
sented in [117]. 

The matrix element for DIS electro-production of J/ip 
mesons has been calculated in [118], keeping the full Q 2 
dependence as well as the full polarization state of the J / ip 
meson. For studying J/ip meson polarization properties 
we calculate the p_i_,v>- and Q 2 -dependences of the spin 
alignment parameter a [19, 61]: 



a(u>) 



da/ckv — 3 da l / du 
da/duj + dah/dbj 



(29) 



where a (ctl) is the production cross section for (longitu- 
dinally polarized) J/ip mesons and with or Q 2 sub- 
stituting ui. The parameter a controls the angular distri- 
bution for leptons in the decay J/ip — > n + (in the J/ip 
meson rest frame): 



dr(j/ip 



) 



dcosO 



1 + a cos 2 9. 



(30) 



Fig. 13 shows the spin alignment parameter a(p±^) cal- 
culated in the region 0.4 < z < 0.9 (a) and 0.4 < z < 1 
(b), with z = E^p/E 7 in the p-rest frame, in comparison 
with experimental data taken by the ZEUS collaboration 
at HERA [119]. Curves 1 corresponds to calculations in the 
collinear approach at leading order using the GRV gluon 
density, curves 2 corresponds to the k± -factorization cal- 
culations with the BFKL JB [120] unintegrated gluon 
distribution with Q 2 ~ 1 GeV 2 . 

We note that it is impossible to make definite con- 
clusions about the fej_-factorization approach considering 
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Fig. 13. The spin alignment parameter a(p±^), which is cal- 
culated in the region 0.4 < z < 0.9 (a,c) and in the region 
0.4 < 2 < 1 (b,d) at y/s = 314 GeV, m c = 1.4 GeV and 
^Iqcd = 250 MeV. Curves 1 corresponds to calculations in the 
collinear approach at leading order with the GRV gluon den- 
sity, curves 2 corresponds to the k± -factorization calculations 
with the JB umntegrated gluon distribution. 



the polarized J ftp production cross section because of the 
large uncertainties in the experimental data. The large 
additional contribution from the initial longitudinal po- 
larization of virtual photons weakens the effect of ini- 



tial gluon off-shellness even more. However at low Q 2 < GeV 2 , 50 GeV < W < 225 GeV, 0.3 < z < 0.9, p* T \ > 1 



Fig. 14. The spin alignment parameter a(p 2 ± ^) at fixed values 
of Q 2 for 40 GeV < W < 180 GeV, z > 0.2, M x > 10 GeV at 
Vs = 314 GeV, m c = 1.4 GeV and /Iqcd = 250 MeV. Curves 
1 and 2 correspond to the calculations as in Fig. 12. 



The CO matrix elements responsible for the non- 
perturbative transitions in eq.(31) are related to the fic- 
titious CO wave functions, that are used in calculations 
in place of the ordinary CS wave functions. In Fig. 15 we 
present a comparison between our theoretical calculations 
and experimental data collected by the HI collaboration 
at HERA [121] in the kinematic range 2 GeV 2 < Q 2 < 100 

* 2 



1 GeV such contributions arc negligible. This fact should 
result in observable spin effects of final J ftp mesons. As 
an example, we have performed calculations for the spin 
alignment parameter a as a function of p]_ ^ at fixed values 
of Q 2 for 40 GeV < W < 180 GeV, z > 0.2, M x > 10 GeV 
at HERA. 

The results of our calculations at fixed values of Q 2 = 
0.1, 1, 5, 10 GeV 2 are shown in Fig. 14. Curves 1 corre- 
spond to calculations in the collinear approach at leading 
order using the GRV gluon density and curves 2 corre- 
spond to the fc^-factorization calculations with the JB 
unintegrated at Qq = 1 GeV 2 . We observe large differ- 
ences between predictions of the leading order of the color 
singlet model with the GRV gluon density and the k±- 
factorization approach at low Q 2 < 1 GeV 2 (Fig. 14). 

Therefore more accurate measurements of polarization 
properties of the J/ tp mesons will be an interesting test of 
the ^-factorization approach. 

For the production of J/ip particles in the framework 
of the CS model the relevant partonic subprocess is 



7 + 5^ 3 Si[l]+3. 



(31) 



When the CO qq states are allowed, there appear addi- 
tional contributions from the following partonic subpro- 
cesses: 

7 + .g^ 1 5o[8]+ 5 , 3 Si[8]+.g, 3 Pj[8]+g. (32) 



GeV 2 . 

The effect of the different evolution equations (BFKL 
JB or DGLAP "derivative of GRV" , for a detailed de- 
scription see [1]) which govern the evolution of gluon den- 
sities is found to be as large as a factor of 2 in the produc- 
tion cross section. This is illustrated by a comparison of 
dash-dotted and dashed histograms in Fig. 15. A similar 
effect is connected with the renormalization scale /i 2 in the 
running coupling constant a s (^L 2 ). The calculations made 
with ii 2 = fc 2 and ^ 2 = to 2 ± are represented by the dash- 
dotted and dotted histograms in Fig. 15. Note that the 
setting fif = k\ is only possible in the fcj_-factorization 
approach. In this case, a s (k±) was fixed at a s =l if the 
formal running value was greater than 1, and it was set to 
zero if k± < Aqcd- The contributions from the 2 — > 1 CO 
subprocesses are cut away by the experimental restriction 
z < 0.9. Turning to the 2^2 CO contributions, one has 
to take care about the infrared instability of the relevant 
matrix elements. In order to restrict the 2 — > 2 subpro- 
cesses to the perturbative domain, we introduce the reg- 
ularization parameter q 2 e g. The numerical results shown 

in Fig. 15 are obtained with setting q^eg = 1 GeV 2 and 
m c — = 1.55 GeV. 

The results shown in Fig. 15 are obtained with the 
non-perturbative CO matrix elements of Ref. [103]. If the 
values extracted from the analysis [117] were used instead, 
the contribution from the CO states would be a factor 
of 5 lower. One can see that, irrespective to the partic- 
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Fig. 15. A comparison between the theoretical predictions and 
experimental data [121] for inelastic J/tp production. Dash- 
dotted histogram, the CS contribution with JB gluon density 
and a a {k]_); dashed histogram, the same with "derivative of 
GRV" and ct s {k 2 L ); dotted histogram, the same with JB gluon 
density and a s (m^ ± ); solid histogram, the sum of the CS and 
CO contributions, with JB gluon density, a a (k 2 L ). 



ular choice of the non-perturbative matrix elements, the 
production of J/ip mesons at the HERA is reasonably 
described within the color-singlet production mechanism 
(with kj_ factorization) and the color-octet contributions 
are not needed. 

4 Selected topics on the current experimental 
status 

In spite of the fact that QCD has been extremely success- 
ful in describing the physics at Q 2 ^> Aqcd, the total 
cross section in deep inelastic scattering (DIS) is domi- 
nated by soft and semi-hard processes which can not be 
described by perturbative QCD. It is thus of fundamental 
importance to provide experimental measurements which 
may give hints to how these processes can be described 
within the QCD framework. 

One of the long standing questions in high energy colli- 
sions is whether significant deviations from the successful 
DGLAP [8-11] evolution equations can be observed at the 
HERA and/or Tevatron colliders. A fundamental question 
is where DGLAP evolution breaks down and emissions, 




An=ln(s / 1) Ar|=ln(x jet /x) 



Fig. 16. Kinematics of hard emissions for di-jets with large 
rapidity separation in hadron-hadron collisions (left) and for 
forward jets/particles at HERA. The maximal measurable jet 
separation at the Tevatron is about six rapidity units. 



not ordered in virtuality, play a significant role. In deep in- 
clastic scattering processes (DIS) at low values of Bjorken 
x it is assumed that the struck parton results from a long 
cascade of parton branchings. Similarly, in hadron-hadron 
collision processes where two jets are separated by a large 
rapidity interval Aq one expects a long partonic cascade 
between them (see Fig. 16). 8 At sufficiently low values 
of x (high values Arf) the DGLAP approximation should 
fail while the BFKL and CCFM approximations should 
be applicable. 

Calculations of inclusive quantities like the structure 
function F 2 (x, Q 2 ) at HERA, performed in NLO DGLAP, 
are in very good agreement with the measurements [83, 
177]. 9 However the interplay of non-perturbative (input 
starting distribution) and the perturbative (NLO DGLAP 
evolution) elements in this calculation makes it impossible 
to decide if parton cascades with strongly ordered trans- 
verse momenta are the dominant mechanism leading to 
scaling violations. 

When exclusive quantities are investigated, the agree- 
ment between NLO coefficient functions convoluted with 
NLO DGLAP parton densities and the data is less satis- 
factory, and for some processes the DGLAP based theory 
fails completely. One example is the cross section of for- 
ward going jets at HERA, which will be discussed below. 
The question therefore is, for which observables the next 
order in the perturbative expansion is enough, and for 
which a resummation to all orders is needed. 

The forward jet production cross section at small 
Bjorken x at HERA and the cross section for jet produc- 

8 Large rapidity interval events in hadron-hadron collisions 
correspond to a subclass of DIS events at small x characterized 
by presence of the forward jet with Xj et = Pjet /p P roton ^> x and 
transverse momentum pr ~ Q. For such events rapidity inter- 
val between forward jet and struck quark Arj ~ ln(xj e t/x) ^> 1 

9 The longitudinal structure function Fl(x, Q 2 ) at small val- 
ues of Q 2 is an exception, see section 4.1.1 
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collincar factori; 

direct 


lation 

resolved 


k±- 

factorization 




LO+PS 


higher order 
NLO (dijet) 


LO+PS 


higher order 
NLO (dijet) 


LO+PS 


HERA obscrvablcs 












DIS D* production 
photoprod. of D* 


ok 


ok [122, 123J 
ok [125,126] 


? 

ok [125] 


? 

no [125] 


ok [123, 124) 

ok [6,127-129,124] 


DIS B production (visible) 
DIS B production (total) 
photoprod. ot B (visible) 
photoprod. of B (total) 


ok [130J 
no [131] 

OK [lOZ, ±OOj 

no [136,133] 


ok [130) 

ok [131] 
? 

no [136,133,126] 


— 
? 


— 
? 


ok [130] 
no [131] 
ok [134, 135] 
ok [134, 135] 


high Q di-jets 

low Q di-jets (cross sec.) 

low Q 2 di-jets (azim.corr.) 

photoprod. of di-jets 


? 

no [139] 
? 


ok [137,138] 
ok )139J 
no [139] 

NLO 3-jet no [139] 
ok [141] 


? 

ok [139] 
? 


? 

no [137, 138, 140J 
? 

no [140, 142] 
ok [141] 


? 

ok [139] 

? 


particle spectra 
energy flow 


no [143, 144) 

no [143, 146, 147] 




ok [145] 

ok [147,145] 




ok [6] 


BJill'tA. small-ic obscrvablcs 












DIS forward jet production 
DIS forward tt production 


no [148-152) 
no [155,156] 


no [151-153J 

? 


ok [149-152,154) 
ok [155,156] 


ok [152, 153) 
? 


ok [6[ 
1/2 [156] 


DIS J ftp prod, 
photoprod. of J f ip 
Jftft polarization 


? 

no [158] 


ok [159] 


? 


? 

ok [160] 

low. stat. [160] 


ok [157,118] 

ok [161,162,118] 

low. stat. [163,118] 


Tcvatron obscrvablcs 


direct 


heavy quark excitation 




D meson prod. 






no 


? 


ok [164,165] 


J/tp prod. 
Xc prod. 

J/tp polarization 


ok [100,102,103,166,167] 
ok [102,103] 
low. stat. [169] 


no 


no 


no 


ok [113,114,168] 
ok [116,168] 
ok [114,168] 


high-p_i_ B prod. 
bb (azim.corr.) 


no [170) 


ok [171) 
ok [171] 


ok [170) 


? 


ok 
ok 


135,164,172-174[ 
174] 


T prod. 


ok [102,103[ 








ok 


168] 


high-p_|_ jets at large Zi?7* 


no 


? 




? 


? 



Table 1. Summary of the ability of the collinear and k ^-factorization approaches to reproduce the current measurements of some 
observables: ok means a satisfactory description; 1/2 means a not perfect but also a not too bad description, or in part of the 
phase space an acceptable description; no means that the description is bad; and ? means that no thorough comparison has been 
made. The label NLO-dijet means the calculation was performed in next-to-leading order for a dijet configuration available for 
example in the DISENT, MEPJET or DISASTER++ programs. LO + PS is a short for LO Matrix Element + parton shower 
calculations as implemented in LEPTO [175] or RAPGAP [176] in the collinear approach and LDC [2-5] or CASCADE [6, 
7] in the k ^-factorization approach. 



tion with large rapidity separation in high energy hadron- 
hadron collisions (Tevatron) have since long been adver- 
tised as an ideal test of the perturbative dynamics [178- 
180,53,70]. More refined theoretical and phenomenologi- 
cal analyses have shown that these tests are not so decisive 
and straightforward, however they remain in the center 
of experimental activity of small x physics. Several mea- 
surements of highly energetic jets (xj e t = Ej et /E p 3> x) 
at large pseudo-rapidities 10 rji a b with transverse energies 
squared Ej, - et of the order of Q 2 have been made by 
both the HI [149,151] and the ZEUS collaborations [150, 
152]. This kind of measurement, originally proposed by 
Mueller and Navclct (so called Mueller- Navelet forward 
jets [178,179,53]) is designed such that DGLAP evolution 
in transverse momentum space is suppressed (Exjet ~ 
Q) while BFKL evolution in .x-space (xj e t 3> x) is en- 
hanced. The measurements showed large discrepancies to 
DGLAP NLO-dijet calculations at low values of Bjorken- 
x, which was taken as an indication of the breakdown 
of the DGLAP approximation and the onset of BFKL ef- 

1(1 The pseudorapidity r\i a b is defined as rji a b = — lntan(#/2), 
with polar angle 9 being measured with respect to positive z- 
axis, which is given by the proton beam (or forward) direction 



fects. However, in the NLO-dijet calculations the scale un- 
certainties are very large. It was also shown [154] that a 
good description of the data can be obtained by consider- 
ing the partonic structure of the virtual photon, which is 
expected to be important for Ef, - et > Q 2 and Q 2 not too 
large. In this approach we are back to the classic DGLAP 
approximation, with simultaneous evolution in transverse 
momentum space from both the photon and the proton 
sides towards the hard scatter and this approach should 
be relevant if chains with at most one local maximum in 
tranverse momentum dominate. HI [151,156] also mea- 
sured single forward particles (tt°) as opposed to jets, al- 
lowing the study of angles closer to the proton direction 
and smaller Bjorken-rr, at the price of a stronger depen- 
dence on the fragmentation process. Studies of the inclu- 
sive cross section for the production of forward particles 
led essentially to the same conclusions as the study of for- 
ward jets. 

Studies of the transverse energy flow provide a com- 
plementary means of investigating QCD processes. Com- 
pared with jet and leading particle studies, measurements 
of transverse energy flow are sensitive to parton emis- 
sions of lower transverse momentum and to the model- 
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ing of both the perturbative QCD evolution and the soft 
hadronization process. In Ref. [156] both types of mea- 
surement (forward particle and energy flow) have been 
merged. The transverse energy flow measured for events 
with forward ir° reveals the range over which the trans- 
verse momentum of the forward parton is compensated. 
We may expect that different models of parton evolution 
lead to different radiation patterns. 

At first look the measurements of the Mueller-Navelet 
jets and jet azimuthal de-correlation at the Tevatron 
should be more promising as a test of non-DGLAP dynam- 
ics compared to HERA due to higher energy and there- 
fore larger phase space open to gluon emissions. While 
at HERA, the separation between the struck quark and 
forward jet can reach up to about four units of rapidity, 
the measurable jet separation at the Tevatron is up to six 
units. Although the BFKL calculations are expected to be 
more reliable at high energies, it should be kept in mind 
that energy-momentum conservation is not fulfilled. Thus, 
effects due to the conservation of energy and momentum 
(consistency constraint) will be significant. 

In Tab. 1 we present a collection of various experimen- 
tal results from HERA and the Tevatron which relate to 
low x physics and we state the result of the comparison of 
these data with NLO DGLAP theory and BFKL and/or 
CCFM evolution schemes. The aim of this section is to 
review some of the items in Tab. 1 in more detail. In the 
next two subsections we review measurements which can 
be well described by the DGLAP approximation and then 
discuss measurements which were especially designed to 
extract a BFKL signal i.e. where the DGLAP evolution 
was suppressed by experimental cuts. 



4.1 Where NLO DGLAP (almost) works 

Several programs exist for the numerical NLO calculation 
of jet observables at the parton level in collinear factor- 
ization. They are known to give comparable results. All 
of them calculate the direct photon contributions to the 
cross sections. Only the JetVip program [181,182] pro- 
vides the additional possibility to calculate a cross section 
consisting of both direct and resolved photon contribu- 
tions, however in the DIS case conceptual difficulties are 
encountered [183] which lead to ambiguous results. 

All the parton level calculations in the collinear ap- 
proach presented in this paper were performed using the 
DISENT program [184] and the CTEQ6M [185] set of par- 
ton distribution functions. The renormalization scale was 
set to [i r = J^P^i an d the factorization scale [if = Q. 11 

There has been some confusion in the literature con- 
cerning the concept of NLO. Formally, whether a calcu- 
lation is leading or next-to-leading depends on the ob- 
servable. LO is then the lowest order in as in which the 
observable obtains a non-zero value, and NLO is one order 
higher in a s . However, sometimes it is difficult to define 
the order in a s appropriate for a specific measurement. 

11 For technical reasons DISENT allows only \Xf = Q or fif — 
const. 



Therefore, in this paper we state clearly to which order in 
a s a process has been calculated. In Fig. 17 representa- 
tive diagrams are shown for NLO calculations of single-jet 
(F2), di-jet and 3-jet processes. It is obvious from Fig. 17 
that NLO dijet calculations do not include all diagrams 
necessary for a NNLO single-jet calculation, as indicated 
in the right column of Fig. 17. 

The calculations summarized above use integrated par- 
ton distributions, convoluted with LO or NLO coefficient 
functions. As is usual in the standard formulation of fac- 
torization, the coefficient functions are on-shcll partonic 
cross sections with subtractions of the singular collinear 
regi ons. The parton distributions are typically in the 
MS scheme, where the partons are integrated over all 
transverse momentum with the resulting ultra-violet di- 
vergences being canceled by renormalization in the MS 
scheme. There is an evident mismatch of approximated 
and exact parton kinematics in such calculations. For the 
infra-red-safe jet cross sections that are the domain of va- 
lidity of the calculations, the factorization theorem ensures 
that the calculations are consistent and valid. 

The transverse momenta of the partons entering the 
hard scattering can be seen as being generated by two 
sources: the intrinsic (primordial) transverse momentum, 
which reflects the Fermi motion of the partons in the 
hadron, typically of the order of one GeV, and the 
transverse momentum generated by the QCD evolution 
(DGLAP or BFKL/CCFM/LDC), which can reach large 
values, even in DGLAP up to the factorization scale. 
Therefore, for more exclusive components of the cross sec- 
tion, it is better to use suitably defined unintegrated dis- 
tributions and off-shell parton kinematics. For one treat- 
ment along these lines that is specifically designed to treat 
NLO corrections in the context of showering Monte-Carlo 
event generators, see the paper of Collins and Zu [186]. 



4.1.1 The longitudinal structure function Fl(x,Q 2 ) 

The longitudinal structure function Fl(x,Q 2 ) is domi- 
nated by the gluon density at large enough Q 2 . In the limit 
of small Q 2 there is no phase space for strong ordering in 
virtuality (DGLAP will not work) and unphysically nega- 
tive values for Fl(x, Q 2 ) are obtained in some calculations 
for Q 2 < 1 GeV 2 [187]. However in the fcj_-factorization 
approach there is no strong ordering in virtuality and 
therefore the parton evolution may generate arbitrarily 
small fc^(down to an artificial cutoff) which means that 
the parameterization is valid over the full range in k± . In 
Fig. 18 the structure function Fl(x,Q 2 ) as measured by 
HI and ZEUS [188] is compared to calculations using k±- 
factorization in the framework of Ref. [189]. The uninte- 
grated gluon density was taken from CCFM (J2003 set 1). 
Shown for comparison is another unintegrated gluon den- 
sity obtained from the derivative of the integrated gluon 
density (here GRV [99] is used) and the prediction for 
Fl obtained in the collinear DGLAP approach using the 
MRST2002 [190] parton densities. 
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Fig. 17. Schematic picture of the diagrams which are included in NLO calculations for single-jet (F2 ), di-jet and 3-jet processes. 
Also indicated are example of diagrams which are not included (since of higher order in a 3 or do not contribute to the n-jet 
observable) 
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Fig. 18. T/ie longitudinal structure function Fl(x, Q 2 ) as mea- 
sured by HI and ZEUS [188] compared to different calculations, 
including also the collinear DGLAP approach. 
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4.1.2 Single inclusive jets at HERA 

Jets have been studied extensively at HERA and other col- 
liders. These measurements have shown that at sufficiently 
large transverse momenta and/or momentum transfers the 
NLO QCD theory based on the DGLAP approximation is 
in excellent agreement with the data. To judge how well 
this approximation works let us mention that the deter- 
mination of the strong coupling constant a s from recent 
HI [191] and ZEUS [192] jet measurements are not only 
in perfect agreement with the world average value but are 
also in precision comparable to LEP measurements. An- 
other example of a jet measurement fully compatible with 
NLO theory in the collinear approach is the measurement 



Fig. 19. Inclusive jet cross section daj^t/dEr in different 
ranges of r/iat integrated over the region 5 < Q 2 < 100 GeV 2 
and 0.2 < y < 0.6. The data are compared to NLO dijet calcu- 
lation (DISENT) and to the predictions from Cascade. 



of dijet angular distributions [193] performed by the DO 
Collaboration. The result of the data-theory comparison is 
an exclusion limit on the quark substructure which is com- 
petitive with many LEP results. In spite of this spectac- 
ular success of the QCD theory in the collinear approach 



20 



Jeppe Andersen et al.: Small x Phenomenology - Summary and status 



> 

2 

2 10 2[ 

H 
LU 

■a 

1 10 
b~ 

■a 



5 < Q < 10 GeV 2 



10 



- • H1 Data 



10 



-2 



10 



-3 



10 



10 < CT < 20 GeV' 



_ 20 < Q < 35 GeV 2 _ 35 < Q < 70 GeV' 



NLO dijet./i,= (,' 
NLO dijet /i,=£ 
LO dijet 
Cascade set 1 
Cascade set 2 



10 



10 



10 



70 < Q < 1 00 GeV 



10 



E T (GeV) 



Fig. 20. Inclusive jet cross section dajet/dEr for the forward region 1.5 < T)i a b < 2.8 in different ranges of Q 2 . The NLO-dijet 
(DISENT) and the LO dijet calculation as well as the prediction Cascade are shown. 



one should keep in mind that there are regions of phase 
space, where the description of the data is less than sat- 
isfactory. It is the aim of this subsection to localize these 
regions and observe patterns characteristic for a possible 
failure to describe the data by NLO theory in the collinear 
approximation. 

In a recent HI paper [194], NLO calculations of the 
inclusive jet cross sections, using the DISENT program 
[184], were confronted with high statistics data. The kine- 
matic range considered in this analysis was constrained by 
the conditions 5 < Q 2 < 100 GeV 2 and 0.2 < y < 0.6 (the 
latter condition leads to a reduction of photoproduction 
background). Jets are defined using the inclusive fcj_ clus- 
ter algorithm [195,196] in the Brcit frame 12 and selected 
by the requirement Exjet > 5 GeV. 

Fig. 19 shows the inclusive jet cross section as a func- 
tion of the transverse jet energy Exjet m different re- 
gions of the pseudorapidity r\\ a b'- in the backward region 
— 1 < Vlab < 0.5, the central region 0.5 < r\i a b < 1.5 and 
the forward region 1.5 < r]i a b < 2.8. The measured cross 
sections, which extend over four orders of magnitude, are 
compared to the calculations obtained in the collinear and 
fcj_-factorization approaches, respectively. 

While there is good agreement between the data and 
the NLO-dijet calculation in the backward region for all 
Er,jet values, discrepancies are observed for more for- 



The Breit frame is defined by 2xp + q = 0, where x is 
Bjorken scaling variable, and p and q are the proton and the 
virtual photon momenta, respectively 



ward jets with low Ex. jet- In the lowest Ex, jet range (5 < 
Ex, jet < 20 GeV, for r]i a b > 1.5), the assumed renormal- 
ization scale uncertainty (Ex,jet/2 < Mr < 2Ex,j e t) does 
not cover the large difference between the data and the 
calculation. In Fig. 20 the forward region from Fig. 19 is 
studied in bins of Q 2 , showing that discrepancies to NLO- 
dijet calculations are most significant at small Q 2 and 
small Ex. jet values. The factorization scale uncertainty is 
estimated by changing /J/ = y^Q 2 to jif = 8.4 GeV, which 
is the average jet transverse momentum. The correlation 
of large NLO /LO corrections and high sensitivity to renor- 
malization scale variations with poor agreement between 
data and QCD predictions strongly suggests that the in- 
clusion of higher order (e.g. NNLO or resolved photon 
component) terms in the QCD calculations is necessary 
in order to improve the description of the data. 

The predictions obtained in the fc^-factorization ap- 
proach, supplemented with the CCFM unintegrated gluon 
densities, as implemented in Cascade, are in reasonable 
agreement with the data. Especially the forward region 
(Fig. 20) is reasonably well described. The quality of the 
various approaches to describe the data can be seen in 
Tab. 2, where we quote the x 2 /nci/, both for the NLO- 
dijet and the Cascade calculations. 

It is interesting to quote in the above context the re- 
cent ZEUS measurement on inclusive jets [198] presented 
in Fig. 21. The cross section for jets reconstructed in the 
laboratory frame with the inclusive k± algorithm is com- 
pared to a calculation in NLO (here 0(a s )). When going 
from small towards large values of r]j et the description of 
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NLO-dijct 
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RAPGAP 










J2003 


dir 


dir + res 


parton density 
factorization scale [i 2 


CTEQ6M 
70 GeV 2 


CTEQ6M 

Q 2 


set 1 


set2 


CTEQ6M 


CTEQ6 
Q 2 + P± 


M + SaS 

q 2 + 4p1 


da/dEt (in bins of »j) (cf. Fig. 19) 

da/dE t (in bins of Q 2 for 1.5 < )) < 2.8) (cf. Fig. 20) 


12.8 


13.2 


25.5 


4.0 


23.7 


1.3 


8.6 


3.9 


13.6 


17.3 


6.0 


13.2 


2.1 


13.2 


da/dAr] (cf. Fig. 23) 


40.1 


40.9 


116.8 


37.7 


66.9 


22.6 


46.7 


S = 1n 2 _U^*,*,Q 2 )^* (cf - Flg ' 24) 


17.8 


15.7 


23.2 


3.9 


3.3 


2.6 


1.7 


forward jets HI p± > 3.5 GcV (cf. Fig. 25(a)) 


8.9 


17.0 


2.7 


4.7 


10.8 


4.4 


0.3 


forward jets HI p± > 5 GeV (cf. Fig. 25(6)) 


5.7 


11.2 


1.9 


2.3 


6.7 


2.6 


0.7 


forward jets HI prel p ± > 3.5 GeV (cf. Fig. 26) 


1.7 


6.4 


1.3 


1.1 


4.4 


0.9 


1.7 


forward jets ZEUS p± > 5 GcV (cf. Fig. 27) 


28.9 


38.4 


19.2 


9.5 


27.1 


20.1 


16.8 



Table 2. Comparison of^/ndf obtained from comparing different predictions to the data. For the NLO-dijet calculation with 
the DISENT program the renormalization scale was set to fj, r = ^2 k± , the CTEQ6M [185] and SaS [197] parton distribution 
functions of the proton and photon, respectively, are used. 



the data by the NLO calculation becomes worse. The rea- 
son for this is that the a° contribution goes to zero and 
the 0(a s ) calculation becomes essentially the LO contri- 
bution, since for a fixed Q 2 and x (or y), rjj e t is fixed in 

an a ° calculation, simply given by ^ In . The range 

used in the measurement is Q 2 > 25 GeV 2 , y > 0.04 
which, for the lowest Q 2 gives a maximum r\j et of about 
0.8. For larger Q 2 the maximum rjj et is a bit larger, but the 
suppression of the a® contribution is still visible in Fig. 21 
around r\j e t = 1. Beyond this, the lowest order contribu- 
tion is dominated by a* and the " NLO" calculation in the 
figure becomes leading order. As can be seen from Fig. 21, 
even DGLAP type Monte Carlo models (here Lepto) give 
a rather reasonable description, if further parton radia- 
tion is included via parton showers. Thus this comparison 
shows the need for higher order corrections but not nec- 
essarily a need for any BFKL contribution. 



4.1.3 Inclusive di-jets at HERA 

The measurement of di-jet production, which is less in- 
clusive compared to the measurement described before, 
might provide a stronger test of the NLO QCD calcu- 
lations in the collinear factorization approach, as it in- 
volves more observables. Experimentally, possible devia- 
tions from the DGLAP approach can best be observed by 
selecting events in a phase space regime, where the main 
assumption, the strong ordering in fc_Lof the exchanged 
parton cascade (Fig. 22), is no longer strictly fulfilled. 
This is the case at small x. The parton configurations not 
included in DGLAP-based calculations might contribute 
significantly to the cross section. Moreover, with respect to 
the photon-proton center-of-mass system (hems), the two 
partons produced in a hard scattering process (Fig. 22) 
are no longer balanced in transverse momentum. Events 
coming from calculations beyond 0(a s ) will lead to a sit- 
uation where the two hard jets are no longer back-to-back. 
The excess of such events is expected to be higher for a 
BFKL scenario compared to DGLAP, due to the possibil- 
ity of hard emissions in the parton evolution provided by 
the non-ordering in k±. 

Di-jet production in deep inelastic ep scattering was 
investigated in the region of low x (10~ 4 < x < 10~ 2 ) 
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section da/dr/j e t for the inclusive phase space compared to ARI- 
ADNE (CDM), LEPTO(MEPS) and NLO Disent m order 
0{a\) Lower part: Relative difference of the measured inclu- 
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and low Q 2 (5 < Q 2 < 100 GeV 2 ) [139]. Jets were recon- 
structed in the hems using the fc^-algorithm. The min- 
imum transverse jet energy of 5 GeV was required 
and an additional requirement on the most energetic jet 
^Tmax > 7 GeV (in the hems) was added to avoid a 
scenario for which NLO-dijet predictions become unreli- 
able [199,200]. In Fig. 23 the triple differential inclusive 
di-jet cross section in bins of Bjorken-a; and Q 2 as a func- 
tion of the distance | Ai]* | between jets is presented. The 
data are compared to NLO-dijet predictions. The NLO- 
dijct calculation with /i 2 = Q 2 falls well below the data. A 
better description over the full phase, including the regime 
of very low x, is obtained using /i 2 = 70 GeV 2 . It should 
be noted however, that even with (i 2 f = 70 GeV 2 the the- 
oretical uncertainty due to scale dependence of the NLO 
calculation is rather large so again no strong statement 
about the DGLAP approximation for dijet production at 
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Fig. 22. Generic leading order diagrams for di-jet produc- 
tion in ep scattering: the variables kj_i, Xi and Zi denote the 
transverse momenta, the longitudinal energy fractions and the 
fractional energy in the splitting, respectively, and p±i are the 
transverse momentum of the radiated gluons. 



low x can be made on the basis of the cross section mea- 
surement alone. In Fig. 23 the data are also compared 
with the predictions using the fc^-factorization approach 
in Cascade. The quality of the data description is again 
quoted as x 2 /ndf for both approaches in Tab. 2. 

Further insight into small- a; dynamics may be gained 
from inclusive di-jet data by studying the behavior of 
events with a small separation in the azimuthal angle, 
Ac))* , of the two hardest jets as measured in the hems, 
as proposed in [201-204]. Partons entering the hard scat- 
tering process with negligible transverse momenta, as as- 
sumed in the DGLAP formalism, lead mainly to back-to- 
back configurations of the two outgoing jets with A(j)* = n. 
Higher order QCD processes lead to azimuthal jet sep- 
arations different from w, however, the effect might be 
smaller than in the case of the BFKL and CCFM evolu- 
tion schemes. In the above quoted di-jet analysis [139] the 
jet azimuthal de-correlation was studied using a variable 
which has been proposed by Szczurek et al. [204] : 



S = 



J*N 2 ^ et (Ar,x,Q 2 )dAr 



Jo 80 N 2 _ jet (A<P*,x,Q 2 )dAr 



-,0<a< 180 (33) 



with a being a parameter for the A<f>* distribution. Its ad- 
vantage in comparison with the direct A<jf measurement 
(see e.g. analysis of the Tevatron data, subsection 4.2.3) 
is its better stability against migrations. For the data pre- 
sented in Fig. 24 a = 120 was chosen. This choice is 
mainly dictated by the limited jet energy resolution which 
may result in an incorrect choice of the two most energetic 
jets. Fig. 24 presents the ^-distribution as a function of 
x in bins of Q 2 . The measured value of S is of the order 
of 5% and increases with decreasing x. This rise is most 
prominent in the lowest Q 2 bin. On the contrary, the NLO 
-dijet QCD calculation predicts 5-values of order 1%, sev- 
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Fig. 23. The triple differential inclusive di-jet cross section in 
bins of Bjorken-x and Q 2 as a function of the distance \An*\ 
between the di-jets compared to NLO-dijet calculation (DIS- 
ENT, solid line) and predictions from Cascade fdas/ied and 
dotted line). 



eral standard deviations below the data, and show no rise 
toward small x. Here the NLO calculations are performed 
in the on-shell limit (see discussion in section 4.1), ne- 
glecting the transverse momentum coming from the QCD 
evolution. Therefore only the 0(a 2 s ) part of the matrix 
elements gives a significant contribution for A(f>* ^ 180 . 
The calculation of NLO-3jet is in much better agreement 
with the data (shown in [139]) which is NLO for the S 
variable, but it still fails to describe the rise towards small 
x. Since Monte Carlo generators, like RAPGAP, include 
the effects of the finite transverse momentum of the incom- 
ing partons via parton showers, it is not surprising, that 
they come much closer to the data than the naive NLO 
calculation ignoring the off-shellncss of the incoming par- 
tons. This shows, that care has to be taken by applying 
fixed order partonic calculations to exclusive observables. 

The CCFM evolution as implemented in CASCADE 
[41] describes the data reasonably well (Fig. 24), but this 
is also true if a resolved component of the virtual photon 
is added, provided that a rather large scale fi 2 = Q 2 + 'ip 2 L 
is chosen (to get large enough resolved contribution). 
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Fig. 24. Ratio S of the number of events with a small az- 
imuthal jet separation (a < 120 ) of the two most energetic jets 
with respect to the total number of inclusive di-jet events, given 
as a function of Bjorken-x and Q 2 . The data are compared to 
NLO-dijet calculations (DISENT, solid line) and predictions 
from Cascade. 



To conclude this section let us summarize its main 
points: 

— For the inclusive jet cross section in DIS, the NLO- 
di-jet description starts to fail when jets become more 
and more forward 

— the worsening of the description is accompanied by 
increasing theoretical uncertainty due to scale depen- 
dence, indicating that the NNLO terms may be more 
important in the forward direction 

— NLO-di-jet calculations describe di-jet cross section in 
DIS data very well, down to XBj = 1CP 4 in the central 
region of rapidity if fi 2 = 70 GeV 2 is used. For [i 2 = 
Q 2 the description is much worse. 

— the largest differences between NLO theory for the in- 
clusive jet and di-jet cross section and the data are 
observed in the small- ir, small Q 2 region. 

— the azimuthal jet de-correlation in di-jet DIS data is 
not described by NLO-dijet calculation, which is ef- 
fectively LO for that observable. The NLO-3jet calcu- 
lation is in better agreement with the data, but still 
at small x is not sufficient. The CCFM evolution ap- 
proach is consistent with the data, but so is LO ME 



+ DGLAP parton shower provided that the resolved 
photon contribution is taken into account with a scale 
given by Q 2 + Ap\. 

As in the case of cross sections, the largest discrep- 
ancies for azimuthal de-correlation are found in the 
small-x, small Q 2 region. 



4.2 Where NLO DGLAP does not work 

4.2.1 Forward jets in DIS 

Measurements described as "Mueller forward jets in DIS" 
[178-180] were especially designed to search for non- 
DGLAP evolution signatures. The following conditions 
were required to suppress the DGLAP and enhance the 
BFKL evolution (with x jet — E jet /E p ): 

— a high energetic jet with an energy fraction Xj e t S> xsj 
to enhance BFKL evolution x-space 

— a high enough transverse momentum Exjet of the jet 
to ensure that perturbative calculations are valid e.g. 
E T .jet > 3.5 GeV 

— Exjet ~ Q in order to suppress the DGLAP evolution 

At HERA, the requirement of Xj et /xBj to be large results 
in typical jet angles of a few degrees with respect to the 
forward (proton) direction. Due to the unavoidable beam- 
pipe hole in the detector, the acceptance is limited to jets 
with an angle larger than, for example, 7° in the HI detec- 
tor. At smaller angles the jets are insufficiently contained 
in the detector and the experimental separation from pro- 
ton remnant fragments might be difficult. As the jet ap- 
proaches more and more the forward direction its profile 
gets thicker and more asymmetric and a large fraction 
disappears down the beam hole. In fact, the observation 
of broadening of the jet profile leads the ZEUS collabora- 
tion to restrict the forward jet analysis to pseudorapiditics 
rjjet < 2.6 corresponding to the limiting angle ©j e t > 8.5°. 
The criterion which determines the minimum acceptable 
jet angle is a satisfactory description of the jet profile. 
Obviously, the jet profile and separation of the remnant 
fragmentation depends on the jet algorithm. In the ZEUS 
analysis [150] the cone algorithm was employed. In prin- 
ciple an algorithm like the fc^-cluster algorithm, which is 
not based on geometry, should be less sensitive to detector 
edges, and the separation of remnant fragments should be 
easier in the Breit frame. 

The condition E^j^/Q 2 m 1 is essential to suppress 
the DGLAP evolution in direct photon interactions. Due 
to limited statistics a compromise has to be found. In 
practice an interval is defined around this central value of 
Erjet/Q 2 - The requirement Exjet ~ Q leads to another 
experimental challenge: to reconstruct jets of the smallest 
possibly transverse momentum, forward jets at smallest 
possibly Q 2 (but still in perturbative region) and hence 
smallest possible x are required. 

It should be noted, that at HERA energies the above 
cuts restrict the phase space not only for DGLAP but 
for any type of evolution. At HERA the range between 
the hard scattering and the forward jet covers about 
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HI cuts 


ZEUS cuts 


K 


> 11 GeV 


> 10 GeV 


2/e 


> 0.1 


> 0.1 


^Tjet 


> 3.5 (5) GeV 


> 5 GeV 


Viet 


1.7 - 2.8 


< 2.6 


■E-T jet/ Q 


0.5 - 2 


0.5 - 2 


X jet 


> 0.035 


> 0.036 


Breit 

.Fzjet 




> 


X 


0.0001 - 0.004 


0.00045 - 0.045 



Table 3. Summary of the different selection criteria used to 
define Mueller forward jets in the HI [151] and ZEUS [150] 
measurements. 



4 rapidity units, limiting the number of hard emissions 
to about 3-4. Therefore, it may be that there is not 
enough phase space for a BFKL-DGLAP discrimination. 
In Fig. 16 (right) a typical Feynman diagram for forward 
jets and particles is shown. 

In Tab. 3 we present cuts used in the HI [151] and 
ZEUS [150] forward jet analyses, performed with the cone 
algorithm. In spite of small differences of the selected 
phase space, it is clear that cross sections at E^jet > 

5 GeV are compatible (see Fig. 256 and 27). 



Recently the HI Collaboration performed a new mea- 
surement of the forward jet cross section [149] using much 
higher statistics and applying cuts almost identical to 
those applied in [151] (5 < Q 2 < 75 GeV 2 , E T . jet > 
3.5 GeV, 7° < 6 jet < 20°, x jet = E jet /E p > 0.035, 
0.5 < E^j et /Q 2 < 2). In this analysis the jets were re- 
constructed using the inclusive /c^-algorithm. The cross 
section for forward jet production [151,149] as a function 
of x is shown in Figs. 25, 26 and 27. The measurements are 
up to a factor of two larger than the prediction based on 
0{a s ) (and also C(a 2 s )) QCD calculations in the DGLAP 
approach. Such parton level calculations are compared in 
Figs. 25, 26 and 27 with the measurement. Also shown is 
a comparison with different unintcgratcd gluon densities 
implemented in Cascade [6,7], which shows the sensi- 
tivity of the predicted forward jet cross section on the 
details of the unintegrated gluon density. It is interesting 
to note that also including the non-singular terms in the 
CCFM splitting function (J2003 set 2) leads to reason- 
able agreement with the measurements in Figs. 26 and 27. 
For all other distributions e.g. da/dp± j et the pattern of 
agreement/disagreement is similar. The level of agreement 
with the measured cross section of forward jet production 
is given in Tab. 2. 

The forward jet cross section was also studied using 
the BFKL formalism [12-14]. In particular Kwiecihski, 
Martin and Outhwaite (KMO) in ref. [205] used a mod- 
ified LO BFKL equation, supplemented by a consistency 
constraint which mimics higher orders of the perturbative 
expansion, to describe the inclusive forward jet cross sec- 
tion. The KMO model describes the data well, however 
the predicted cross section is very sensitive to the input 
parameters in particular to the infrared cut-off and the 
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Fig. 25. The cross section for forward jet production [151] at 
the hadron level, as a function of x for (a) Exjet > 3.5 GeV 
and (6) Erjet > 5 GeV. Also shown are the predictions from 
LO and NLO dijet calculations as well as predictions from Cas- 
cade. 



scale of a s . Thus the model has rather large uncertain- 
ties in the normalization of the cross section, whereas the 
shape of the distribution in Bjorken-x is expected to be 
more stable. We will come back to the KMO calculation 
in the next section. 

The choice of the jet algorithm has quite an effect 
on the measured cross section, as we can see comparing 
Fig. 29 (data from [149], k± jet algorithm in Breit frame) 
and Fig. 28 (data from [148], cone jet algorithm in lab- 
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Fig. 26. The cross section for forward jet production[149] at 
the hadron level, as a function of x for Erjet > 3.5 GeV. Also 
shown are the predictions from LO and NLO dijet calculations 
as well as predictions from Cascade. 



Fig. 27. Forward jet cross section [150] as a function of x in 
the kinematic region given in Table 2. Shown are LO and NLO 
dijet calculations together with predictions from Cascade. 



oratory frame). The cross sections come out different at 
hadron level due to the choice of the jet algorithm. 

Including a contribution from resolved virtual photons 
(as done in Rapgap res [176]) leads also to a reasonable 
description of the forward jet data. It should be noted 
however, that the predictions of the model are sensitive 
to the renormalization and factorization scales. The Rap- 
GAP package allows a choice of renormalization and fac- 
torization scale, and in Figs. 28 and 29 the predictions 
are presented for two different choices, /i 2 = Q 2 + p\ and 
H 2 = Q 2 + ^p\i where p± is the transverse momentum 
of the partons taking part in the hard scattering pro- 
cess. The errors (mainly systematic) are large and rea- 
sonable agreement with the data would still be achieved 
for a scale of Q 2 + \p\ . Note, that for a correct descrip- 
tion of the azimuthal de-correlation by Rapgap the same 
large scale has to be employed (see subsection 4.2.2). How- 
ever it seems that the two different forward jet measure- 
ments prefer different choices of the scales: in Fig. 29 
the forward jet data are well described with a renormal- 
ization scale fi 2 . = Q 2 + p\ while the forward jet data of 
Fig. 28 lie between the predictions using [i 2 = Q 2 +p± and 
A* 2 = Q 2 + ^Pj_- Both calculations use Rapgap and the 
same (CTEQ6M and SaSld) proton and photon PDF's. 

Before coming to the end of this subsection, let us 
comment on possibilities of a new type forward jets mea- 
surements, which open with the advent of high statis- 
tics data. Obviously we can go to higher Q 2 and higher 
Xj et so that the ratio Xj et /xBj would remain large. An 



example of such a scenario is Q 2 > 16 GeV 2 ,p± j e t > 
6GeV,x jet > 0.05,0.1 < y < 0.7. The cross section 
calculated using RAPGAP and CASCADE is shown in 
Fig. 30. It is 3-5 times lower than measurement with cuts 
presented in Tab. 3, but is certainly measurable at the 
level of 100 pb _1 and has several advantages. The jets 
with higher p± and higher Xj et are cleaner, we can ex- 
pect smaller systematic errors due to the uncertainty of 
the calorimeter scale and detector corrections. Further- 
more in the region of higher Q 2 the resolved component of 
the photon is suppressed, therefore the ambiguity between 
CASCADE-like and RAPGAP-like descriptions may van- 
ish. Another possibility was considered by Kwiecifiski et 
al. [206] who studied deep inelastic events containing a for- 
ward photon as a probe of small- a: dynamics. The great 
advantage is that such a measurement is no longer depen- 
dent on the hadronization mechanism. At an integrated 
luminosity of around 1 fb _ we can expect about 300 
BFKL-like events within the following phase space cuts: 
20 < Q 2 < 30, k 2 ± > 5 GeV 2 , 6> 7 > 5° where k 1± and 6> 7 
are the transverse momenta and the angle of the forward 
photon, respectively. 13 The DGLAP theory prediction is 
about 3.5 times lower. This process seems to be measur- 



It is necessary to impose an isolation cut on the photon 
to suppress background from 7r°'s produced within outgoing 
quark jet. Experimentally one requires that within isolation 
cone around photon the energy deposit is below few percent of 
the photon energy. 
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Fig. 29. The predictions of the RAPGAP Monte Carlo 
(dir+res) at two different scales compared to the data of [149] 
(k± jet algorithm) . 
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Fig. 30. Forward jet cross section as a function of Bjorken 
x in new cut scenario designed for high statistics data : Q 2 > 
16GeV 2 , pTjet > 5 GeV, x jet > 0.05,0.1 < y < 0.7. 



2 < CT < 4.5 GeV z 
HI preliminary 

DIR + RES, Q 2 + 4p,2 
CCFM (CASCADE) 
DIR 

mod. LO BFKL 



4.5<Q Z < 15GeV" 




15<Q"<70GeV / 



Fig. 31. The cross section for forward ty° production as a func- 
tion of x forp* Tlt a > 2.5 GeV. Also shown are the predictions 
from various Monte Carlo calculations. 



able at HERA 2, however background from 7r°'s may turn 
out to be large. 

4.2.2 Forward 7r° mesons 

HI recently measured single forward ir° meson produc- 
tion [156]. This new measurement triples the number of 



7r°'s in comparison to previously published data [151], 
allowing the measurement of more differential distribu- 
tions and additional final state observables. The analysis 
is restricted to the kinematic range 2 < Q 2 < 70 GeV 2 , 
5° < W < 25°, = En/Ep > 0.01 (lab. system) and 
Pt v > 2-5 GeV (hems). The differential cross section as a 
function of Bjorken- a; for different regions in Q 2 is shown 
in Fig. 31. It should be noted that this measurement covers 
a range in x down to 4 ■ 10~ 5 . The prediction of a DGLAP 
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based Monte Carlo (Rapgap dir) is well below the data, 
whereas a reasonable description is obtained when the re- 
solved virtual photon contribution is added (Rapgap res) . 
It should be also noted, that a rather large factorization 
and renormalization scale Q 2 + Ap\ has to be used in this 
case. Surprisingly, Cascade (all sets, but only JS2001 is 
shown) falls below the data at small x values. The fact 
that Cascade, which provides good description of the 
forward jet production, fails to describe the forward 7r 
production at small x is interesting in itself. It may indi- 
cate that quark initiated cascades and final state cascades 
(gluon splitting into quark pairs), both missing in present 
Cascade generator code, play important role in the for- 
ward ir° production. In Rapgap both processes contribute 
significantly to the forward ir° cross section, influencing 
both the scale for string fragmentation (string invariant 
mass) and the string composition (quark vs gluon frag- 
mentation). The final effect is such that Rapgap is able 
to produce significantly more forward 7r°'s. It is interesting 
to note, that the parton to hadron fragmentation usually 
viewed as a complication of the partonic picture of deep 
inelastic collisions, here may serve also as the indicator 
of the underlying parton dynamics. It should be stressed 
however, that there is no direct contradiction in the data: 
discrepancies in the ir° cross section arise in the region 
of x which is mostly beyond the reach of the forward jet 
measurement. It is interesting to note that the previously 
mentioned BFKL calculation of the forward jet cross sec- 
tion [205] is consistent with the forward ir° mesons mea- 
surement (dashed-dotted curve) . The BFKL prediction for 
the 7T° cross section was obtained by the convolution of the 
parton distribution of KMO [205] with the KKP (Kniehl, 
Kramer, Potter) fragmentation function [207]. 

We expect that different initial state cascade dynam- 
ics should lead to different radiation patterns and there- 
fore to a different transverse energy flow. The trans- 
verse energy flow in the hadronic center of mass sys- 
tem, j^dE^,/d(rj* — r]*), in events containing at least one 
forward ir° is presented in Fig. 32 where 77* gives the 
pseudorapidity of the pion (in the hems). The spectra 
are presented in three intervals of the 7r° pseudorapid- 
ity: — 1.25 < 7]* < —0.25 (closest to proton direction) , 
-0.25 < 77* < 0.25 and 0.25 < rj* < 2.0 (farthest from 
proton direction). 

The QCD-based approaches all describe the transverse 
energy flow around the ir° but give different predictions 
in the current region. The resolved photon picture gives 
a reasonable description of the spectra whilst the CCFM 
approach overestimates the transverse energy flow when 
the forward ir° is closest to the proton direction (top left). 
The direct photon model gives the worst description of 
the data. It predicts a transverse energy flow which rises 
strongly with increasing Arj* and shows a peak at large 
values of pseudorapidity difference. This effect becomes 
less pronounced with increasing pseudorapidity of the for- 
ward 7T°. The differences between the models can be qual- 
itatively understood as a consequence of the ordering cri- 
teria of the parton cascade implemented in various Monte 
Carlo generators. 




Fig. 32. The distribution of transverse energy as a function of 
pseudorapidity difference in different intervals of tv° pseudora- 
pidity. Predictions of three QCD-based models are shown. 



4.2.3 Di-jet production at the Tevatron 

The jet production data at high energy hadron-hadron col- 
liders can also be used to test parton evolution dynamics. 
The production of exactly two jets is described at LO by 
an a 2 calculation as being back-to-back in azimuthal angle 
and having their transverse momenta balanced. Higher or- 
der processes involve the radiation of additional partons, 
which will upset this correlation and additional soft radi- 
ation in higher order processes will decrease the correla- 
tion further, leading to a smearing of the /^-distribution. 
Perturbative QCD has been successful in describing di-jet 
production up to next-to-leading order, whereas higher 
order contributions have to be accounted for by parton 
shower models. Since the production mechanism may in- 
volve more than one hard interaction scale, a different 
treatment of the parton radiation, such as BFKL, might 
be needed. 

The DO experiment has studied events in which two 
jets, widely separated in rapidity, have been identified. 
Due to their uniform calorimetric coverage of ± 4 units 
in rapidity this experiment is well suited for such an in- 
vestigation. Jets were defined using a cone algorithm and 
E t jet > 20 GeV. In a multi jet event the two jets mostly 
separated were chosen for the analyses provided one of 
them had E t > 50 GeV, this in order to avoid any trigger 
inefficiency. 

If (cos(7r — A<j))) is plotted as a function of the rapidity 
separation between the observed jets, then one would ex- 
pect to observe a decrease in this variable as the rapidity 
separation increases, simply because the phase space for 
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additional radiation increases. As shown in Fig. 33 the DO 
experiment [208] also observes a linear decrease with the 
pseudorapidity interval, well described by the HERWIG 
Monte Carlo. The JETRAD Monte Carlo, which pro- 
vides a NLO-dijet calculation, predicts less de-correlation 
at large rapidity gaps. The BFKL calculation by [70], 
valid for large s, on the other hand gives much larger de- 
correlation effects, although we note that in this analysis, 
large effects from the constraint of energy and momentum 
conservation have been ignored in the BFKL evolution. 
In fact, if these are taken into account as describe in Sec- 
tion 3.3 the BFKL prediction is in much better agreement 
with data [55]. 

Recently the DO measurements of Mueller-Navelet jets 
at the Tevatron have been discussed in detail by Andersen 



et al. [59] , therefore we restrict ourselves to quoting their 
main conclusions: 

— Definitions of the momentum fractions used by DO and 
some of the acceptance cuts imposed spoil the correct- 
ness of the procedure to extract the effective BFKL in- 
tercept from the data. Especially the implemented cut 
on the maximum allowed transverse momentum of jets 
invalidates a BFKL analysis based on the asymptotic 
behavior of the BFKL prediction. Such a cut will of 
course always be implicitly implemented by the con- 
straint in energy at a given collider, necessitating a 
detailed analysis as described in section 3.3. 

— As the cuts on the transverse momenta of trigger jets 
were chosen equal, the fixed NLO QCD calculations 
of both the total dijet rates and the azimuthal de- 
correlations are plagued with large logarithms of per- 
turbative, non-BFKL origin. 

The constrained phase space at the Tevatron for dijets 
with large rapidity separation puts severe limits on the 
phase space for mini-jets (contributing to the BFKL evo- 
lution). The phase space constraint prohibits the rise in 
cross section with increasing rapidity separation (simply 
because the decrease in the PDFs is faster than the in- 
crease in the partonic cross section), but other observables, 
like the angular correlation of dijets, still get large BFKL 
corrections. The LHC promises to be very well suited for 
a study of effects from the BFKL evolution. 



4.3 Experimental Conclusions and Outlook 

The measurements of forward jets and particles are sen- 
sitive to the dynamics of parton evolution. Several QCD- 
based approaches have been confronted with the data. It 
has been shown that NLO-dijet DGLAP calculations fall 
well below forward jet data. The forward jet cross sec- 
tion is, however, well described by a DGLAP based Monte 
Carlo which includes a resolved photon component. Simi- 
larly, results obtained using the BFKL and CCFM evolu- 
tion schemes are compatible with the data. The measure- 
ment of the forward ir° cross section leads essentially to 
the same conclusions for the region of Bjorken-a; covered 
by the forward jet measurements. For the lowest values 
of x i.e. those beyond the reach of the forward jet mea- 
surements, the CASCADE Monte Carlo generator fails to 
describe the data, possibly due to missing quark initiated 
cascades. The comparison of results from various models 
seems to indicate some sensitivity to the fragmentation 
method used to connect the parton and hadron levels. 
Study of the transverse energy flow associated with for- 
ward 7r°'s seems to favor the DGLAP direct + resolved 
approach. 

However, the present measurements at HERA were 
mainly restricted by two factors: the available center-of- 
mass energy and the geometrical acceptance of the detec- 
tors, requiring the forward jet to lie between: 2 ^ rj ^ 3. 
The dijet measurements arc described best with a different 
scale (Rapgap) or unintegrated gluon density (Cascade) 
than the forward jet measurements (cf. see Tab. 2). This 
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shows that indeed new effects are seen: if the forward jet 
cross section is extended to a range of r)j et up to 6 units 
(as proposed in the proposals for a continued HERA3 pro- 
gram [209,210]) the difference compared to DGLAP be- 
comes even more significant. 



5 Conclusions 

On the theoretical side, significant progress in understand- 
ing small- a; effects has been made. The soft region (k± ~ 1 
GeV) has been clearly identified to have significant influ- 
ence on hadronic final state observables. With the consis- 
tent treatment of the scale in a s and including the non- 
singular terms into the CCFM splitting function, a nec- 
essary step forward to a serious application of the CCFM 
small- a; evolution equation has been taken. 

The question of gauge invariance of the whole k±- fac- 
torization approach in general and also the question of 
gauge invariance of (integrated or unintegrated) PDFs has 
been clarified further. 

Many new measurements in the area of small- a; physics 
have been made public, and the interest in a better un- 
derstanding of small- a; effects is very clear. New mea- 
surements indicate the need to go even beyond O(af) 
if calculations are performed in the collincar factoriza- 
tion approach. On the other hand, these effects are auto- 
matically included in k±- factorization, which makes k±- 
factorization an important tool for studying higher order 
corrections. It was shown that, irrespective of the partic- 
ular choice of the non-perturbative CO matrix elements, 
the production of J/ ip mesons at HERA can be reasonably 
well described within the color-singlet production mech- 
anism (within fc^-factorization) and color-octet contribu- 
tions are not at all needed. More accurate measurements 
of polarization properties of the J /if) mesons will be an 
interesting test of the fcj_-factorization predictions. 

Still, results from the experiments at HERA and the 
Tevatron have not yet provided unambiguous evidence for 
new small- a: effects. Including higher orders in the calcula- 
tion according to the collinear approach and/or including 
the concept resolved (virtual) photons seems to mimic also 
new small- a; effects. In order to unambiguously identify 
small- a; effects at e.g. HERA, it is necessary to increase 
the angular coverage of the experimental setup towards 
the proton direction, as has been shown in the proposal 
for an extended HERA running beyond 2006, the so-called 
HERA 3 scenario. Since a correct description of the small- 
x dynamics is essential for the understanding of QCD at 
high energies, and also for any asymptotically free field 
theory it is of great importance to continue and extend 
the experimental and theoretical efforts. 
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